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1 .0 Sutima ry 

This design note documents the entry-to-touchdown Space Shuttle 

Orbiter simulation capability that has been incorporated into the 

GEMASS (ref. 9) Subprogram 33 (3-DOF), The simulation makes use of the 

guidance as outlined in the 19 Nov. 1976 SSOFT Level C FSSR 

(ref, 1). A digital autopilot interfaces between GEMASS and the 

guidance-. Vehicle attitude is determined by use of .the ability 

of GEMASS to integrate differential equations in addition to the 

equations of motion. Vehicle aerodynamic characteristics are 

obtained from an OV-102 Aerodynamic Data Tape and control' surface 

deflections required to trim the vehicle and trimmed aerodynamic 

coefficients are determined internally. Several indicators to 

allow evaluation of .subsystem performance have been included; as 

is the ability of the user to activate any of several error 

dispersion sources. The performance of the simulation compares 

well with more sophisticated 6-DOF simul ations’ as will be shown 

later in this' note in a comparison with SSFS (ref. 10). The simulation is 

currently available on the 1108-EXEC2 system and will soon be 

available on the 1110-EXEC8 system. Execution time on EXEC2 

system is approximately 25% that of SSFS on the EXEC8 system 

which makes this a very economical simulation tool. This v/ork 

was conducted under contract Number NAS 9-14960 Task Order B0705. 
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2.0 Introduction 

The entry-tb-touchdown- Space Shuttle Orbiter simulation was 
developed to provide an economical tool for subsystem performance 
and. evaluation studies for the Integrated Entry Task, (EX43) of 
the Engineering Analysis Division (EAD) of NASA. Previous 
simulations used in. this -task were capable of simulating the 
flight of the Onbiter front entry interface (400 K ft) to M = 2.5 
ysing the Analytic Drag Control '(ADC)- closed-loop guidance. An 
open-1 ooD guidance routine (Cl,CT) was also available, but an open- 
loop guidance cannot achieve acceptable touchdov/n. conditions. 

Since actual Shuttle flight wil-1 be control 1 ed by closed-loop 
guidance techniques, extending the open-loop capability 
further was- deemed to be of .little value. Therefore, work was 
u'rtdertakeh tb extend the closed-loop capability of the GEMASS 
simulation. 

Add-ition of the Terminal Area Energy Management (TAEM) guidv-nce 
allowed: extension' of the closed-loop simulation to M = .9. At this 
point the TAEM guidance would, initiate modulation of the speed- 
brake. The trimmed aerodynamic data tables being used by GEMASS 
were created only for the nominal speedbrake deflection schedule 
and adding tables- of trimmed aero data for various speedbrake 
positions or including speedbrake- 'as an independent- variable 
in the tables would have been prohibitive in terms- of runtime, 
core requirements, or program modifications;' Therefore,' it was 

i 

necessary to add the ability to calculate trimmed ^aero data 
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internally to the simulation and access the Aerodynamic Data 
Tape directly. To complete the simulation, the Approach and 
Landing (A/L) guidance was added to provide the commands necessary 
to deliver the vehicle to acceptable touchdown conditions. 

Two existing ■ quasi-degrees of freedom in angle of attack and 
bank angle were used to simulate the transient response of the 
vehicle to the guidance commands. The digital autopilot determines 
the accelerations and rates of these flight variables and the 
rates and' accelerations are integrated by 6EMASS along with the 
3-DOF equations of motion. 

Several additions were made to the simulation to allow evaluations 

of subsystem performance that would not be normally attainable 

in 3-DOF programs. Chief among these is a model of its the Thermal 

Protection Subsystem (TPS) with calculations of surface and 

backface temperatures and required insulation weights and 

thicknesses. Control surface hinge moments computations are 

done by the simulation; in addition, the user has the option of obtaining 

punched output that can be used in the sonic boom program to 

obtain ground overpressure levels for the entry traa'ectory . A payToad 

bay venting pressure indicator was also added. 

The ability to apply several off-nominal dispersion sources has 
been included in the simulation. Among these are increments to 
the aerodynamic coefficients, Guidance and Navigation errors, and 
atmospheric dispersions. In this latter area, two modes of 
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dispersion value application to the neutral atmosphere are 
available. Winds may be applied with or without activating 
an alpha-wind correlation routine. Omni -direction winds. or 
a statistical wind model v/ith gust and shear buildup may be 
used. 
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3.0 DISCUSSION 

The routines that distinguish the entry-to-touchdown simulation 
from the other GEMASS subprograms are those concerned with guid- 
ance, navigation, targeting, and the determination of trimmed 
aerodynamic data. Each of these is discussed in some detail 
below, and all interface with GEMASS via the digital autopilot. 
Obviously, this autopilot (DAP6EM) does more than the function 
implied by its name (this will be discussed later in some detail) 
as it was convenient to perform the requisite additional func- 
tions there. Enclosure 1 shows a block diagram of the information flow 
between GEMASS, the digital autopilot, and the other routines. 

The routines on the right hand side of the page are calledM'n 
order from top-to-bottom of the page. The top line going into 
each routine shows the inputs required; the bottom line coming 

i 

from each routine shows the outputs to DAPGEM. None of the 
routines on the right hand side of the page access each other, 
although some (such as targeting and navigation) furnish inform- 
ation which is used in routines called later, 0n”ly one guidance 
routine is used at a time and previously used guidance routines 
cannot be reentered i-n the same case. DAPGEM computes all input 
variables that are not available in GEMASS; i.e. pitch and roll 
angle. Others (like relative velocity) are accessed by the 

r 

routines from the GEMASS A-array. No flight variables are modi- 
fied by any of the routines; new variables are created (such as 
by , coordi nate transformations) which are required for the execu- 
tion of the routines. 
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Besides serving the functions of providing and accepting informa- 
tion from the routines of the entry-to-touchdown simulation, 

DAP6EM provides GEMASS with aero coefficients in the form 
required, provides the rates and accelerations of angle of attack 
and bank angle which are integrated by GEMASS to simulate two 
additional degrees of freedom, provides additional useful informa- 
tion to the user, and assures that the guidance routines are 
called in the proper sequence and at the proper rate. DAPGEM 
is accessed by GEMASS through the DERIV routine which computes 
the external forces to be used in the integration of the equations 
of motion. For this reason DERIV is within the loop of the inte- 
gration technique being used (Runge-Kutta or Adams-Moulton) and 
the digital autopilot is called each integration cycle. However, 
the guidance, navigation, and targeting routines must be cycled 
at specific rate, DAPGEM and. the user required input discussed 
later in Implementation assure that this is done. Thus, 

DAPGEM regulates accessability to the guidance, navigation, 
and targeting. 

The remainder of this discussion will deal with the routines 
mentioned on Enclosure 1 in greater detail with emphasis on 
user required input and useful output which may be desired 
during execution of the entry-to-touchdown simulation. 

Other routines that are available for use in Subprogram 33 such 
as the TPS model, sonic boom output, open-loop guidance, 
dispersion source modeling and other minor subroutine mods will 
be discussed. 
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DERIV 

This, routine computes the external forces to be used by GEMASS in 
the equations of motion. It requires that the user put an integer 
value into A{1360) to select the type of guidance to be used. A value 
of less than three (3) is required to access the open-loop guidance. See 
the discussion on CLCT for the options available. A value of g.reater 
than two (2) is required to access the closed-loop guidance for the. 
entry-to-touchdown simulation. See the following discussion on DAPGEM 
for the options available. 

Besides the choice of which guidance option to use, two additional 
sets of computations were put into DERIV for the entry-to-touchdown 
simulation. The viscous interaction parameter (See reference 2, Sec. 3.11) 
•is computed here. This parameter is required during the early portion 
of entry as one of the Independent- variables in determining the aero 
data components. A(1351) to A(1354') are us'ed. to compute the viscous 
interaction parameter which is stored in A(1353). 

The trajectory heating rate and load is computed by the equation 
in reference (3). • A(IOIO) to A(1024) are used to store the inter-r 

mediate results. The heat rate and load are stored in A{591) and A(592), 
respectively. These are required for use of the TPS subroutine. 
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DAPGEM 

This routine calls the navigation, targeting, guidance and aerodynamics. 

routines. It sends each the required computed inputs, accepts coinnands 

and interfaces with the mainstream of GEMASS (through the subroutine 

DERIV) to obtain body orientation, control surface deflections, and 

aerodynamic force coefficients. DAPGEM calls the three guidances in 

their normal sequence. A previous guidance cannot be re-entered unless 

a new case is initiated. Furthermore, the start of each guidance must 

be coincident with the start of a new phase. This is because a change 

in integration step size is required for each guidance and this can 

only be accomplished by starting a new phase. (Note that starting a 

new case also starts a new phase.) For the user to start a simulation 

in the entry guidance, A(1360) should be set to 3. In TAEM or A/L 

guidance A(1360) should be 4 or 5 respectively. The simulation will 

execute-until a user defined stop condition is met and will then 
1 

continue with the next guidance phase. This process will continue 
until the case or job is terminated. The section on Impliraentat ‘on 
will innumerate what is required to make the proper phase changes »nd what is 
required at the start of each phase or case. 

The entry-to-touchdown 'simulation does not attempt to model the baseline 
Flight Control System (FCS) because of the core and runtime (due to the 
small time steps required) required. The difficulty in employing 
only those parts of the FCS applicable to a 3-DOF simulation would 
also be prohibitive. To simulate the transient response of the 
vehiclje to the attitude and control surface deflection commands from 
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the guidance routines; second-order control equations are contained in 
DAPGEM. The equations determine the angle of attack and bank angle 
accelerations and rates. These quantities are integrated by GEM ASS 
(along with the 3-DOF equation of motion) to determine the body 
attitude. Thus the simulation has two additional quasi-degrees of 
freedom. The additional equations integrated are: 


angle of attack: 
bank angle: 



a 





The second order control equations are:. 

ci - (a^ - a) - a 

(fi = K, (4 - 4) - K: 4 

where: K^» K,jj, are gains (user selected) 

a^= commanded angle of attack 
= commanded bank angle 

When TAEM- and Approach and Landing guidance are being used, DAPGEM 
converts the normal acceleration commands into alpha commands. 


Maximum rates on alpha and phi can be controlled by the user (A(1216) 
and A(1236), respectively), as can maximum accelerations (A(1217) and 
A(1237). Gains on the rate term in the control equation are controlled 
by A(1214) for alpha and A(1234) for phi. Gains on the angular error 
are input in A(1215) for the alpha error and A(1235) for the phi error. 
The user is required to put the initial angle of attack and bank angle 
into A{1211) and A(1231), respectively. 
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DAPGEM also produces a speedbrake rate (limited by the design opening 
and closing rates of 5 dps and 9 dps) from a first order control 
equation which is integrated by GEMASS to produce a speedbrake deflection. 

The speedbrake control equation has the form: 

^SB "" “ “^SB^ / 

where: = commanded speebrake position 

c 

At = current time step (integration step size) 

The integration equation used by GEMASS to determined the speebrake 
position is 

The user is require to input a speedbrake schedule in a table which 
is followed down to M = .9. After M = .9 the TAEM and A/C guidance 
provide the speedbrake command. 

DAPGEM also retracts the body flap when this flag is set by AUTLND. 

The user can control the angle to which the body flap retracts by 
A(1226). The rate of retraction is 3 dps. When the landing gear 
deployment flag is set, DAPGEM produces a post-multiplier for the 
landing gear Increments to the aerodynamics. The rate of deployment 
can. be controlled by the user with A(T228). 
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EGRT 

This routine computes targeting parameters required by the guidance. 
Inputs necessary are the target longitude and geodetic latitude which 
are put into A-array 106 and 107, respectively, the runv/ay azimuth, 
runway ID and altitude above the Fischer ellipsoid (A(1165), A(n80), 
A(1192)). A redesignation capability is provided. The redesi.gnation 
flag (A(T166)) is s.et and a velocity at which redesignation is to be 
done (ACnsy)) is chosen. The redesignation latitude, longitude, 
altitude, azimuth, and runway ID (the runway' number) are input (A(1161), 
A(1162), A(1163), A(1165), A(1185)). Only one redesignation per 
entry is allowed. 

The approach geometry, in the form of the X coordinate of way point 2 

(the nominal touchdown point), the radius of the heading alignment 

circle, and the distance from the V/P2 to X-coordinate of the HAC 
» 

(A(1169), A(1171), A(1172)) are required. EGRT provides the range 
to WP2 and the azimuth error to the ADC guidance (C0NGID),_ It 
provides the vehicle position and velocity in runway coordinates to 
the TAEM and .Approach and Landing guidance. Those two routines have 
the range and heading computations done internally. 
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GIDNAV 

This routine provides several navigation parameters to the guidance; 

These are altitude above the landing site., altitude rate, drag, and 
total acceleration. (A(1196), A(1198), A(1199), and A(1194) respectively). 
No special user inputs are required. These quantities can be included 
in the output if. desired. 
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CONG ID 


This routine contains the ADC guidance as detailed in the 19 NOV 
1976 FSSR (Reference 1 ). Required user inputs (mission depen- 
dent or MX constants in the FSSR) are the commanded angle of 
attack schedule in A(llOl) to A(1141) and parameters which define 
the drag-velocity profile to be flown in A(1142) to A(1160*) and 
A(1177). Altitude rate feedback inputs are in A(1173) to A(1176). 

An. option in the altitude rate feedback scheme al1ov>;s for considera- 
tion of whether the FCS is the automatic or Control Stick Steer- 
ing (CSS) mode. Even though the entry- to-touchdown simulation 
has no FCS, the capability of simulating the guidance change 
resulting from the FCS mode has been retained. If A(1164) is 
set to zero, the normal (automatic) mode calculations are used; 
if A(1164) is set to 1, the CSS mode calculations are used. The 
nominal velocity for Entry/TAEM -interface is i-nput via a( 1200); 
when the simulation reaches a velocity less than this, A(llOO) 
is set to 4. This "flag" may be tested as a stop condition for 
the ADC guidance phase. The ADC guidance returns the commanded 
bank angle and angle of attack to DAPGEM. 

A brief description of the function of the aforementioned 
A-arrays is given in Enclosure 2. For a more detailed description 
of the ADC guidance the reader should see Section 4.6 of Reference 1 
Table 4. 6. 2. 1-2 of this document also gives a description of the 
MX constants and references subsections that describe their 
function in more detail. Where possible the Fortran names in 
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CONGID correspond to the .variable names in the FSSR document. 
Where this was not possible, the Fortran names have been made as 
much like the variable names as possible and the actual variable 
names are contained in parentheses in Enclosure 2 , This same 
procedure was followed in the TAEM and A/L guidances. Table A-1 
of the FSSR contains the recommended OFT-1 values of the MX 
constants . 
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TAEMG 

This routine contains the 19 NOV 1976 FSSR Teminal Area Energy 
Management guidance for OFT-1. It sends a commanded speedbrake 
deflection (for M < .9), a normal force command and bank angle 
command to DAPGEM. Transition criteria in TAEMG are tested to 
determine when the guidance should be exited. When an exit 
condition is met A(llOO) is set to 5; this may be used as a stop 
condition for the TAEM guidance phase. Mission dependent (MX) 
constants that define the altitude and dynamic pressure profiles 
and the energy boundaries are now input by the user via A(1801) 
to A(1850). A brief description of thes« A-arrays. is given in 
Enclosure 2. For a more detailed descriptron of the TAEM gui- 
dance the reader should see Section 4.7 of Reference 1.- Table 
4. 7. 2. 1-3 of this document also gives a description of the MX 
constants and references subsections that describe their function 
in more detail. Table B-1 of the FSSR contains the recommended 
OFT-1 values of the MX constants. 

The quantities EMEPCl, EMEPC2, ENCl , and ENC2 are 2x2 matrices 
that should be loaded column-wise. The rows define values to be 
used with each steep glide slope (GAMSGS) and the column value 
used is determined by the range- to-go. The first steep glide 
slope is used for a vehicle weight of less than 250K Tbs and the 
second is for greater than 250K lbs. All other variables (except 
XA) are vectors of dimension two where the value used is deter- 
mined in the same manner as the steep glide slope. The quantity 
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XA'is a vector of dimension two. that defines the aim points on the 
runway. Initially, the first aim point is used but the user may 
switch aim points, at any time based on one of several conditions 
being met. These conditions are input in A(1851) to A(1854). 

They are (in order); 1) time from start of phase greater than a 
preset value, 2) velocity less than a preselected value, 3) alti- 
tude above the runway less than a preselected value, 4) range-to- 
go less than a preselected value. No values are built in, so the 
user should input sufficiently large or small values to prevent 
the aim point change if not desired, i.e., values that cannot 
possibly be attained during the simulation. The aim point change 
can occur only once per simuTation. The flag for the aim point 
change is A(1855). Zero means no change has occurred; one means 
the change has occurred. This fl.ag is carried into the A/L 
guidance. The quantity TGGS in the NOV FS3R need not be input. 

It is the tangent of the steep glide slope angle and is computed 
internally to -lessen the user input load. 

Several quantities computed by the TAEM guidance have been put 
into the A-array (A(1242) to A(1249)) for potential output. These 
are defined in Table B-2 of Reference 2'. 
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AUTLND 

This routine contains the 19 Nov. 1976 FSSR Approach and Landing 
guidance. Mission dependent (MX) constants that define the reference 
trajectory to be flown are novA input by the user via A(1871) to A{1892). 
A brief description of the function of these quantities is given in 
Enclosure 2. For a more detailed description of the A/L guidance the 
reader should see Section 4.8 of reference 1. Table 4. 8. 2-3 of this 
document also gives a description of the MX constants and references 
subsections that describe their function in more detail. Table C-1 
of the FSSR contains the recommended OFT-1 values of the MX constants. 
The steep glide slope angle and the runway aim point are the same as 
those used in the TAEM guidance and need not be input again. The 
remarks about the aim point change made in the TAEMG discussion apply 
here as well. If an aim point change v/as made in TAEM, it can not 
and need not be done again here. The quantities DSBCTD, SIGMA, and. 

VREF are vectors of dimension tv/o where the first value is associated 
with the steep glide slope for a vehicle weight of less than 250 K 
lbs and the second is for greater than 250 K lbs. The quantities HK, 

XK, R, and XEXP are 2X2 matrices where the rows define the values 
to be used with each steep 'glide slope and the columns define values 
associated with the aim point change flag. The first column is used 
if no change has occured; the second is used if the change has occured. 
These matrices should be loaded columnwise. 
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Two additional user inputs are available for simulating the change 
in the guidance effected by operating the FCS in the automatic or the 
manual mode. If A(1201) and A(1202) are set to zero the guidance 
computations for use of the automatic pitch and roll control by the 
FCS are done. If they are 1, guidance computations simulating .manual 
pitch and roll control are done. 

Several quantities computed, by the A/L guidance have been pu.t 
into the A-array (A(1246), A(1247), and A(1278) to A(1280)) for 
potential output. For a brief explanation of these .^.and other output 
quantities see Enclosure 2. For further explanation the FSSR should 
be consulted. AUTLND' sends the same commands as TAEMG to DAPGEM. 

In addition, it sets flag commands that cause the body ‘flap to begin 
to retract and the landing gear to begin deployment (A(1225) and 
A(1224) respectively). A touchdown flag (A(1232)) is set when the 
altitude of the wheels above the runway is zero. 
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WIND 

This routine performs the alpha -wind correlation documented in 
ref. 4, For details of how the actual angle of attack is computed, 
this reference should be consulted. A{ 121 2) is the navigated angle 
of attack and is used by the control equations in DAPGEM and the' 
guidance routines. A(1310) is computed by WIND and is always the 
actual angle of attack outside of the guidance, with or without 
winds. This allows for no impact in other areas of GEMASS (i.e. 
sonic boom, TPS, aero tables, etc.) where actual angle of attack is' 
required. Note that A(51) controls whether or not this subroutine 
is used. 

The velocity at which to stop the alpha-wind correlation is input 
in A(1207). 
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TAERO 


6EMASS did not have the ability to internally calculate trimmed aero- 
dynamic data which was required to acheive the degree- of fidelity' 
desired of the simulation. Trimmed aero data tables were generated 
off-line as functions of angle of attack> Mach number, and body flap 
deflection for a scheduled speedbrake deflection. Sets of tables 
were prepared far several vehicle center of gravity locations and 
for viscous 'and inviseid aerodynamics. Addition of more Independent 
variables to the tables was not possible since GEMASS is limited to 
three independent variables per table. To add tables to account 
for additional variables would have been prohibitive in terms of core 
storage and runtime. Therefore, it was necessary to add to GEMASS 
the capability of determining the contribution of each component of 
the aerodynamics and determining the trimmed aerodynamics internally 

for the fl'ight conditions desired. This has the additional advantage 
* 

of using the Aerodynamic Data Tape as the source for the tables used • 
instead of card input tables or tapes of card image of input tab’'es 
as was previously done. The files on the Aero Data Tape are set up-'" 
to read by the GEMASS table reading routines. The computation of 
trimmed aerodynamics internally allows extension of the GEMASS 
simulation to touchdown. It allows use of other components of the 
aerodynamics such as speedbrake, landing gear and ground effects 
increments which would be difficult to account for with the previous 
aero data input method. 



This routine obtains all the required aerodynamic components and 
determines the conditions necessary to trim the vehicle (CmTOTAL ' 

It does this through modulation of the body flap and/or elevon. An 
eleven schedule and an initial body flap .deflection (A(1665)) are 
required inputs along with a dynamic pressure (A(1168)) at which to 
start the trim calculations. The process by which the trim elevon 
and body flap deflections are determined is shown in Enclosure 3.. 
Basically, it involves finding a body flap deflection that trims the 
vehicle with the scheduled elevon deflection. If this cannot b.e done 
the body flap is deflected to whichever of its limits (-11.7 or 22.5 
deg) is closest to Cfj, = 0 and the trim elevon deflection is' determined. 
After the elevon deviates from its scheduled value, TAERO attempts to 
return it or move it nearer if possible. Once it- does return, modula- 
tion of the body flap begins again. TAERO will perform the trim 
calculations for and X-Z center of gravity combination '(A(1271) and 
A(1272))» respectively) where X-e.g. is measured in inches aft of the 
forwardmost (Body Sta.t1on 1076.7 or 65% Lg) c.g. position and Z-e.g. 
is measured positive in inches above WL375, TAERO returns trimmed 
and Cq and the required control surface deflections to DAPGEM. 

The dynamic pressure at which the trim calculations start' should be 
2 psf. If an entire entry trajectory is under considera.tion, the 
switch from the viscous to the inviscid aero components should be 
done when the viscous interaction parameter reaches a value of ,005. 
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The alternate table technique normally available in GEMASS may be used 
to transfer to the inviscid tables and A{1358) is reserved for this 
purpose. Similarly, a transfer from the high to low Mach number aero 
table should be done when Mach number reaches a value of 1,25. A{1277) 

Is reserved for this purpose. A(1227) should be used as a post-multiplier 
on the landing gear components. 
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ROUTINES ADDED TO GEMASS 

Several routines have been added to the entry-to-touchdown 
simulation to enhance its usability in performance and subsystem 
verification. 

TPS - This routine contains the aerothermal , weights, and 

backface temperature models as defined in the RI/Houston 
Simplified Weight Synthesis Program (ref 3) through 
revision 4. of that program. A brief description of 
these models is contained in sec 4.8 of reference 5. 
Required user inputs are a flag (A(1400) equals one) 
to access the routine, the initial bare top temperatures 
(A(1434 to 1438) and A(1440 to 1444), all recommended 
to be set to 680), the TPS stop conditions (A(1458 to 
1461)) and the number of panels (A(1651), 25 for 
current models). The first two stop condition inputs 
are A-array locations of flight variables to be tested 
against; i.e., stagnation heating rate (A(591)) and 
altitude (A(108)). The next two are the values these 
flight variables are to obtain before the TPS calculations 
are to b'e stopped. If the stop is desired’ when the 
flight variables are decreasing in value, the values 
should be input v/ith a negative sign. 

The user should also input table control cards 101 to 
106 (respectively) to obtain panel design heat load 
(A{1667)), normal shock viscosity {A(1657}), panel 
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heat ratio (A(1455)), panel transition Reynolds number 
(A(1452)) , panel thickness (A(1656)), and panel area 
(A(1652)), Examples of these control cards are contained 
in Enclosure 5. They access the tabular data tape 
described in Operation and Implementation. Results of 
the TPS calculations are stored in A(140l} to A(1749) 
and all are available for output. A summary page ‘is 
given in the simulation when a stop condition fs reachred 
that gives maximum surface, temperatures , their tiroes of 
occurrance, backfaca temperatures, panel thicknesses, 
weights, and heat loads, and totals. 

SBOOM - The routine can be accessed to obtain punched output 
that can be used in the Sonic Boom program to predict 
overpressure levels. Parameters output are Maoh number, 
altitude, and time rate of change of the variables 
velocity. Azimuth, and flight path angle. Also output 
are flight path angle, azimuth, latitude, longitude, 
bank angle and angle of attack. The routine may be 
accessed by setting A(1950) to 1. The frequency of 
the output is control by A{1951) to A(1998) which are 
input by the user. The first four locations of this 
group are Mach numbers that define where the output 
described above is to start and stop. It is divided 
into two phases. Output can only occur between the Mach 
numbers defined by A(1951) and AC1952) or A(1953) and 
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A{1954). Within each of these phases the output frequency 
can be controlled by the remaining locations, A(1955) is 
the A-array location of the variable that is to be 
tested against for the output check. The remaining 
locations are the values of this quantity at which the 
output should be done. They should be input in the order 
in which they will occur in the simuTation » i.-e., 
ascending or descending.. Tf' threy .are "encountered in 
descending order they should be input v/ith a negative 
sign. 

VECWIN - The vector wind profile detailed in ref (6) has been 
programmed into GEMASS. The A-arrays used for input 
are A(1282) to A(1293). See Enclosure 2 or ref (6) 
for a brief description of the quantities required. 
Associated routines also included are TLKUP and 6ETDAT 
and a common block WND. VECWIIi’s use requires that a 
special Input tape be mounted but table control cards 
are not required since the GEMASS table look-up routines are 
not used. 
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CLCT 

This routine contains three options for open-loop guidance. The 
initial angle of attack is input in A{1310) and the angle of attack 
command A(12T0) is determined from Table 42. A phugoid damping 
term is applied to changes in the command to obtain the actual alpha 
('A(1310)). The three options available to determine the bank angle 
are controlled by A (13,60) : 

A(1300) == 0 ; bank angle is obtained -for a constant q/g trajectory. 
A(1360) = 1 ; constant bank angle from A(1380) 

A(1360) = 2 ; bank angle from Table 41 

To use the constant q/g trajectory, an initial bank angle must be 
input. This angle is maintained until pull-out (y = 0) occurs. 

The heating rate at this time (from reference ll)is 



^ C|_ s Cos <j> 


2 

v/here: q is- in BTU/ft sec 

W is the vehicle weight.' 
is the bank angle, 
is the coefficient of lift. 

K is a constant of proportionality dependent upon 
the entry conditions 

2 

S is the reference area - 2690 ft . 


- 27 - 



Bank angle is then modulated to maintain this heating rate until the 
load factor limit (A(1363)) which is determined from Table 46 is 
reached. This load factor schedule is then maintained for the remainder 
of the trajectory by modulating bank angle. Several A-arrays have been 
reserved for gains in the open-loop guidance. These are A(1304), 
A(1305), A(1312), A(1323), and A(1327), See Enclosure (2) for the 
terms to which these gains are applied. 

A(1361) is input 

by the user for the direction of the bank angle: 0 is for right, 

1 is for l^eft. After the heading is changed from the initial heading 
by 90 deg. the side force coefficient (Cy) is set to zero, but the 
bank angle is maintained. Two options exist for determination of 
the aerodynamic components. If A(1359) is zero, and 6 .^^ are 

determined from Tables 20, 16, and 45 respectively for the inviscid 
data. Tables 47, 48, 49 are reserved for the viscous data. The 
transfer to the inviscid tables should be done when A(1353) reaches a 
value of .005. A body flap schedule in Table 45 must be supplied also. 
If A(1359) is one, the trimmed aerodynamics are determined from- TAERO; 
thus, all the requisite tables must be set up. 
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GEMASS ROUTINE MODS- 


Several routines in the, GEMASS program were modified to be 
compatible with the entry-to-touchdown simulation. Usually these 
require no additional user input but where applicable these will be 
noted. The routines and the extent of their mods are: 

TABLE - As noted in the discussion -on TAERQ., the trimmed data 
tabl es were read into 'cof-e from -a tape. The user had 
to input a number which indicated the C.G. to be used 
and TABLE would select the proper files to be used from 
the tape. This was removed from TABLE since C.G. 
position is taken into account in TAERO. The user 
input is no longer required. 

0UT2 *■ Only 2000 A-arrays are now .printed .out .since the others 
are being used for .data storage. ’This was done to 
keep the A~array as small as possible so the simulation 
would fit in the available core. Note that CT (52) 
should now be set to 2000 and table storage should 
start at 2Q00. 

INTDEP - This routine has the values of the additional 

independent eouations to be integrated. These are 
angle of attack and -bank angle and their rates and 
the speedbrake position. The rates are now stored 
here for phase changes . 
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OUTl 


ATM0S3 


GRAV 

TERMEX 


This routine was modified to allow more output of 
results per printed page. The total number of lines 
possible per page is now 57 including header and two 
blank lines before each block of outout. Normal 6EMASS 
input is required. 

Computation of atmospheric conditions during early 
entry v/as modified to be consistent with pther 
■simulations. 

The call to VECWIN was included in GRAV 
This routine is outside the integration loop of GEMASS 
and is called at the end of every time step. The calls to 
TPS and SBOOM are made here. The hinge moment calculations 
are done here. Required inputs are A{1080) (flag value 
of one required to perform -ca.l^culTat'lons^ ^A(108S;)., 

A(1089), A(1095) (area-arm multipliers for body flap, 
eleven, and speedbrike. See Enclosure 5 or reference 2 
for recommended values) and A(1090). This last Input 
is an increment which is added and subtracted from 
the actual eleven deflection to find the maximum and 
minimum hinge moments that can be anticipated on- the 
elevens due to aileron deflection or off-nominal 
eleven positions. Hinge moments are calculated for 
both the inboard and outboard eleven surfaces. The 
user is also required to provide table control cards 
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to access the fi.les from aero data tape. Enclosure 5 
shows these {Tables 120 to 127). See Enclosure 2 
(A(1080) to A(1096)) for an explanation of quantities 
used in determining hinge moments. A venting pressure 
Indicator (A{1098)) and the local vent pressure 
coefficient (A(1099)) as defined by ref. (7) have 
also been added. 
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4,0 OPERATION AND IMPLEMENTATION 


The GEMASS entry-to-touchdown simulation is contained on each of the 
first two files of tapes X07475 and X06938. A listing *of a sample run 
deck with control cards is shown in Enclosure 5. This run deck is 
sufficient to simulate the OFT-l entry trajectory. The job is divided 
into three phases, one for each guidance mode. The guidance inputs 
required for each phase are shown in the input block for that phase. 

It is not necessary to do this since there Is no overlap in the A-array 
locations of the guidance inputs. Similiarly, a-jcomplete set of the 
output quantities are given for each phase even though this means that 
some of the quantities are repeated unnecessarily. This was done to - 
allow the user to see all the output quantities in each phase. The 
plot input blocks for each phase have been changed to display meaningful 
parameters. The table control cards are input in the first phase and 
not changed.’ 

After execution begins, the first card GEMASS reads sets the position 
in the A’-array where table storage is to begin and the number of 
locations to allocate to table storage. The size of the A-array in 
the entry-to-touchdown simulation is 16000. Since the first 2000 
locations are used for input, output and internal operations the 
remaining 14000 are available for table storage. This is sufficient 
for storage of all aero, TPS, hinge moment, dispersions (see Appendix I), 
schedules, and other user- input tables in their current forms and 
allovs's approximately 1000 locations for future storage requirements. 

The next block of Input is the CT-array. 
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Shown next are the user input tables. In this case, these are the 
eleven and speedbrake schedules (respectively) which are user required 
input. They are input as functions (X -argument) of Mach number. 
Following the tables are the case control, hollerith, and phase control 
cards. In the A-array input block that is next; notice that several 
options (i.e., TPS, hinge moments, etc) have been activated in this 
example. The input block has been divided into sections (as far as 
this is possible) that are required for GEMASS ■'ex’e’cution and those 
that are peculiar to the entry-to -touchdown simulation. 

The comments to the right of each A-array assignment shouTd be 
sufficient to define its purpose. There are a few assignments in this 
block that are of special interest. Note that A(1360) Is set to- 3 to 
start with the ADC guidance. The integration time step (A(53)) is 
1.92 sec. as this is the rate at which the ADjC guidance should be 
sequenced. The number of addftionaT differential .equations to he 
integrated (A(48)) is 5, lie., angle of attack and bank angle rate and 
acceleration and speedbrake rate. Also, the stop is when the guidance 
mode indicator (A(llOO')) becomes 4. The TAEM guidance cannot be 
entered until this occurs. If the user wishes to terminate the ADC 
guidance phase at some point other than that defined by VTAEM, then 
A(llOO) would have to be set to 4 in the phase in which the TAEM 
guidance is to be used. 

After the plots are set up, the table control cards define the source 
and location of the input data tables to fae stored in the A-array. 
Enclosure 6 gives the file locations, size, and arguments of all tables 
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that are currently available for use with the entry-to-touchdown 
simulation. The first page of this enclosure is for a special data 
tape (X14728) which contains many frequently used tables. The 
remainder of the enclosure is for -the special aero data tape (X17990) 
that was created in order to 'run the-simulation. For more information 
on this tape see reference 8. 

The second .phaae (-TAEM) of the simulation also has special input 
requirements. The integration time step must now be set to .96 sec. 
in order to properly sequence the TAEM guidance. The stop condition 
should be when the TAEM guidance transition criteria set A(llOO) to 5. 
A(nOO) must equal 5 for initiation of the A/L guidance. Also, the 
heat rate calculation (A(590)) should be turned off as it will cause 
the simulation to terminate due to an argument out of range for velocities 
less then 1200 fps. It is not required for the TPS routine sines TPS 
has stopped by this time I 

The third phase requires a time step of .16 sec. for proper sequencing 
of the A/L guidance. The stop condition for this phase is touchdown, 
i.e., when A(1232) equals 1. Another possible stop condition could be 
when the vehicle is 18 ft. above the runway {height of the vehicle 
c.'g, above the v/heels). In this case, the runway altitude is 2220 ft. 
above the Fischer ellipsoid so the c.g. altitude at touchdown is 2238 ft. 

A new case could now be initiated starting again with the ADC guidance. 
However, the input constants for the TAEM and A/L guidances must be 
input again. The amount of input required in these phases would be 
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less if the input constants that do not overlap (i.e., do not have to 
be changed from one phase to the next) were all input in the first 
phase of the first case. 

Enclosure ? shows a sample Form 588 and plot request. This provides 
the operator with sufficient information to mount the correct tapes 
and start a GEMASS entry-to-touchdown simulation. If sonic boom 
output .is to be obtained the box for punched output should be checked. 
Runtimes are for complete entry-to-touchdown simulations of one case 
per job. Additional cases would take approximately 6 min. Individual 
phase execution times are; 

ADC guidance - 3.5 rain. 

TAEM guidance ~ 1.0 min. 

A/L guidance ~ 1.5 min. 

An additional 2 min. is required for reading in the program and input. 
The maximum' runtime reflects these requirements plus some contingency. 
The p:ge limit is for printing output from all time steps of a complete 
one case job. The plot number is applicable to the sample run deck of 
Enclosure 5. 
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5.0 RESULTS 

A comparison of a GEMASS entry-to-touchdown simulation of OFT-1 
with a similar trajectory obtained using SSFS (6-DOF) is shown in 
Enclosure 8. Points from the GEMASS trajectory are plotted on the 
SSFS Cal comp plots. The first eleven pages are for the complete 
trajectory; the remainder show in detail the flight variables from 
Entry/TAEM interface to touchd.awn. 

The match between independent flight variables that are related 
to the vehicles position and velocity vector; i.e., dynamic pressure 
and altitude are generally good. Where differences do occur, they 
can be attributed to the two additional quasi-degrees of freedom 
incorporated into GEMASS. While the angle of attach and bank angle 
histories are similar to that of the 6-DOF program, the differing 
response characteristics of the enxry-to-touchdown -stimulation, do 
result in small trajectory differences-. An example of this can be 
seen in the 4000 ft, difference in di vnrange at Entry-TAEM interface 
(Enclosure 8 page 15). Although alpha and bank angle histories were 
similar up to this point, they differed- eno.ugh from the 6-DOF 
simulation to cause this. The downrange error caused the differences 
in bank angle and alpha during the TAEM guidance phase (Enclosure 8- 
pages 20 and 21 ) . 

Finally, it should be noted that the method of determining the 
control surface positions differs between the two simulations. In 
the 6 DOF simulation the attitude commands going to the FCS result in 
control surface deflection commands coming from the FCS and going to 
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the actuators models , The moments which result from this modeling 

determine the vehicle attitude and dynamics. In the 3-DOF 

simulation, the control surface deflections are those required to trim 

the vehicle at its current angle of attack. The eleven history in 

the 3-DOF simulation is the input elevon schedule; the body flap is 

used to trim the vehicle until it begins to retract to its landing 

position. Limits on the elevon and body flap rates are not 

considered. Tn the 6-DOF simulation, the elevon is allowed to 

'deviate up to 1 deg. from its schedule at which time movement of 

the body flap begins to return the elevon to within this 1 deg. 

limit- Rate limits on the surfaces are considered in the 6-DOF 

simulation. Differences in the speedbrake position in the two 

simulations can be attributed to actuator dynamics in SSFS and 

differences in the dynamic pressure betv/een the simulations which 
( 

was caused by angle of attack differences. 
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6.0 - CONCLUSIONS AND FUTURE WORK 

The entry-to-touchdown simuTation capability is operational on 
the Uni vac 1108 EXEC Z system. Although it does not contain an 
Integrated Flight Control System and is a 3-DOF program, its simulation 
capability compares well with more sophisticated 6-DOF programs such 
as SSFS. Program input and usage follow the general GEMASS guide- 
lines and requirements although there are now more options available 
to the user than were formerly. The user has ..available two types of 
attitude commanded guidances and one of these may be exercised in 
any or all of the three major modes as currently defined for the 
baseline Space Shuttle guidance. Simulation of the performance of 
several subsystems may also be done. The program is economical to 
use in terms of runtime which makes it a tool of considerable value 
in doing performance and trade studies. This facet of the program 
becomes particularly evident when multiple case -jobs are submitted. 

Enclosure 9 shows a review of the tapes and the files on those 
tapes that may be used with the entry-to-touchdown simulation. 

Enclosure 10 gives the names and phone numbers of" those programmers 
who should be contacted if there are questions on the simulation. 

The entry-to-touchdown simulation will undergo continual 
modification as new models are added and existing models are updated. 
This will be especially true of the guidance which will undoubtly 
have to be revised as new guidance FSSR's are issued. It may become 
necessary to maintain versions of the simulation for both OFT and 
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Operational Space Shuttle guidance configurations. A tape will he 
created which will contain on separate files all the inputs necessary 
to run any of several trajectories. All the user need to do is 
assign this tape in place of the scratch C-tape and start his job as 
though It were a second case. This process is described in Section 
3.3 of the user's manual. It may be advantageous to add a Flight 
Control System that will more accurately si'mulcite "th‘e “real worT 
control surface deflections and their rates. The simulation will 
soon be operational on the EXEC 8 system to preserve the program 
when and if the EXEC 2 system is finally phased out. It will also 
allow flexibility in usage as each system sometimes experiences 
hardware, software, and turnaround problems. 

As the simulation currently exists on the Univac 1108 EXEC 2 
system there is only limited opportunity for future growth. Less 
than 100 actual core locations are not used for .program or storage. 
The growth limitation may be alleviated by dumming out routines 
that are not required by the user for his simulation. The routines 
to v/hich this may be done, the resulting core mode available, and 
an example of the input necessary to achieve this are shown in 
Enclosure 11. This process may be temporary or new tapes may be 
created with these dummy routines. Note that this last option 
would permanently lose some of the 6EMASS program capability. This 
is especially true of the EXECT 1, EXECT 5, and EXECT 8 routines 
which control several of the other subprograms. This problem could 
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also be solved by further segmentation of the subprogram 33 
segment. This would involve -loading only one guidance routine at 
a’ time. The open or closed loop guidance would not be loaded if 
it were not to be used. This would require remapping the program. 
A disadvantage of this will be an increase in runtime (how much is 
unknown at this time) as the program tape must be read again for 
each guidance; for a multiple case job this would have to be done 
several times. 

As changes are made to the GEMASS program they will be 
reported as updates or appendices to this design note. 
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ENCLOSURE 1 
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Enclosure' (1) 
Page 2 of 4 


A/I 


FLAG 


BF 


FLAG 


"Z 

C6 

D 

h 

h 


LS 


"RW 

L-/D 

LG 


FLAG 


MX 


ADC 


MX. 

MX 

M 

N. 


TAEM 

A/L 


^ZC 


q 

R 


Definition of Variables 

Approach and Landing initiation flag 

Body flap retract flag 

Drag coefficient 

Lift coefficient 

Side force coefficient 

Vertical lift coefficient 

Center of gravity location 
Drag acceleration 
Geodetic altitude 
Altitude above landing site 

Runway altitude 

Lift-to-drag ratio 
Landing gear deployment flag 

ADC guidance mission dependent constants 

TAEM guidance mission dependent constants 
A/L guidance mission dependent constants 

Mach number 
Normal load factor 

Load factor command 

dynamic pressure 

Magnitude of radius vecxor rrom earxn cencer 
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Enclosure' (1) 
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R 

^RW 

RT6 


T/A geom 

h 

Vl 


W 

X, 


a 

a 


c 




a 

AA 

5 


BF 


'EL 


'ELS 


'LG 


Definition of Variables 
Altitude rate 

Position vector in runway coordinates 
Range-to-go 

Target and approach geometry 
Airspeed 

Relative velocity vector 
Inertial velocity vector 
Velocity vector is runway coordinates 
Wind vector 

Viscous intruction parameter 
Weight 

Position vector 

angle of attack 
Commanded angle of attack 

Navigated angle of attack 

Angle of attack acceleration 
Heading error 
Body flap position 

Elevon position 

Scheduled elevon position 

Landing gear deployment ratio 
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Enclosure- (1) 
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Definition of Variables 


^SB 

Speedbrake position 

^SB 

Speedbrake rate 


Speedbrake command 

0 

Euler pitch angle 

X 

Longitude 


Geodetic latitude 

h 

'Bank angle 


Commanded bank angle 


Bank angle rate 


Euler roll angle 
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A-ARRAY FOR ENTRY-TO-TOUCHDOWN SIMULATION 


A-Array 

1000-1005 

1006 

1007-1009 
1010.-1024 
. 1025 
1026-1028 
1029 

1030-1033 

1034 

1035 

1036 

1037 

1038 

1039 

1040 

1041 

1042 

1043-1050 

1051 

1052 

1053 
1054' 

1055 


Definition 

Not used 

Density scale height- ft 
Not used 

Used to compute heat rate and load 
Not used 

X-Y-Z wind component post-multipliers 
Density increment post-multiplier 
Post-multipliers for aero data uncertainties 
Post-multiplier for altitude rate error 
Post-multiplier for range-to-go error 
Not used 

Pressure factor post-multiplier 
Sound speed factor post-multiplier 
Viscosity factor post-multiplier 
Flag for altitude error: 0-off, 1-on 

Altitude error-ft 
Altitude error post-multiplier 
Not used 

Velocity error- fps 

Bank angle error-deg 

Angle of attack error-deg 

Altitude rate error selector: 1-WTR, 2-ETR 

Altitude rate error-fps 
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A-Array 

1056 

1057 

1058-1059 

1060 

1061 

1062 

1063 

1064 - 

1065 

1066-1079 

1080 

1081 

1082 

1083 

1084 

1085 

1086 

1087 

1088 

1089 

1090 


Definition 

Range-to-go error selector: 1-WTR, 2-ETR 

Range-to-go error-NMI 
Not used 

Aero uncertainty series selector: 1-10 

Cq uncertainty 
C[_ uncertainty 

Density factor selector: 1-70°N Hot, 2-ETR Hot, 

3-70°N Cold, 4-ETR Cold 

N/S wind profile selector: 1-N, 2-S 

E/W wind profile selector: 1-E, 2-W 

Not used 

Hinge moment computation indicator: O-Off, 1-On 

Body flap hinge moment-in-lbs 
Inboard eleven hinge -moment-in-Tbs 
Outboard eleven hinge moment-in-lbs 

Inboard eleven hinge mument with positive def. increment- 
in- lbs 

Outboard elevon hinge moment with .positive def. increment 
in- lbs 

Inboard elevon h-i nge. moment , with .negative def. increment- 
in-1 bs 

Outboard elevon hinge moment with negative def. increment 
in-lbs 

2 

Body flap area-arm-ft in 
2 

Elevon area-arm-ft in 

Delta deflection to be considered on the e1evons-deg 
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A-Array 

1091 

1092 

1093 

1094 

1095 

1096 

1097 

1098 

1099 
‘1100 

1101-1110 

1111-1120 

1121-1130 

1131-1140 

1141 

1142 

1143 

1144 

1145 

1146 

1147 

1148 


Definition 

Body flap hinge moment coefficient 

Inboard eleven hinge moment coefficient 

Outboard elevon hinge moment coefficient 

Speedbrake hinge moment coefficient 

2 

Speedbrake area-arm-ft in 
Speedbrake hinge moment-in-lbs 
Not used 

Corrected venting pressure indicator - psi 

Local vent pressure coefficient 

Guidance mode indicator: 3-ADC, 4-TAEM, 5-A/L 

CALPO-constant coefficients of alpha command segments-d'eg 

CALPl -linear coefficients of alpha command segments-°/fps 

CALP 25 quadratic coefficients of alpha command segments- 
/fps^ 

VALP-veloci ty breakpoints between alpha command segments-fps 
NALP-number of alpha command' segments less one 
AK-describes slope betv/een temp, control and eq, glide phases 
AKl-describes slope between drag quadratics 
VA-velocity of zero slope of second drag -quadratic-fps 
VA'1-velocity of intersection of drag quadratics-fps 
VA2-velocity of zero slope of first drag quadratic-fps 
VBl -start velocity of eq. glide phase-fps 
VSl-describes slope of drag profile of eq. glide phase-fps 
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A-Array 

1149 

1150 

1151 

1152 

1153 

1154 

1155 

1156 

1157 

1158 

1159 

1160 

1161 

1162 

1163 

1164 

1165 

1166 

1167 

1168 


Definition 

2 

ALFM-constant~g phase drag~fps 
PREBNK-initial entry bank angle-deg 

EEF4-enesgy value of transition phase drag profile anchor 
point-ftTsec'^ 

DF-drag value of transition phase drag profile anchor 
point- fps 2 

RPT1 -range from anchor point to WP2-NMI 

VTRAN-transition phase start v.e1 oclty-fps 

2 2 

ETRAN-transition phase start energy- ft sec 
VQ-constant of ranging velocity- fps 

651- ro11 smoothing parameter between drag quadratics 

652- ro11. smoothing parameter -between TC&EG .phases 

GS3-ro11 smoothing parameter between EG&C6 phases 

GS4-ro11 smoothing parameter between C6 and transition 
phases 

RTLATD-geod. latitude of redesignation target point-deg 

RTLONG-longitude of redesignation target point-deg 

RTALTD-al titude above Fischer ellipsoid of redesignation 
target point-ft 

(CSS (CSSRYLAMP) - control stick steering flap - 0 automatic 
mode,-l CCS mode 

RRAZ-azimuth of redesignation runway-deg (E of N+, W of M~) 
ITREDS-redesignation flag: 0-no redesignation, 1-redesignate 

VREDS-velocity at which redesignation is to occur-fps 
Dynamic pressure to start trim aero-psf ‘ 
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A- Array 

1169 

1170 

1171 

1172 

1173 

1174 

1175 

1176 

1177 

1178 

1179 
1.180 

1181 

1182 

1183 

1184 

1185 

1186-1188 

1189-1191 

1192 

1193 

1194 

1195 

1196 


definition 

XWP2-X coordinate of way point 2-ft 
RAZ nominal runway azimuth-deg (E of N+, W of N-} 
RTURN-radius of heading alignment circ1e-ft 
DBAR-distance from WP2 to X-coordinate of HAC-ft 
VRDT - velocity to start alt. rate feedback - fps 
ACNl - time constant for alt. rate feedback - sec. 

DDLIM •- max. delta drag, for alt. rate feedback - ft/sec 
ZKl - gain for alt. rate feedback 
D230 “ Initial valve of ref. drag at VBl - ft/sec? 
IPHASE-gui dance phase indicator all modes 

LODV-commanded L/D 

IRWID(RWID) Runway identification index (RW number) 

2 

D23-reference drag at VBl -fps 

2 

DREFP-reference drag accel'erati'on-fpd 
RDTRF-reference h-dot with HSCDDOT correction-fps 
ALDREF-reference L/D 

Redes-lgnation runway identification index 

RG-vehicle X-Y-Z in RW coordinates-ft 

VG-vehicle veloci-ty 'components tn RV/ ,coordinates-fps 

RTE1 -altitude of nominal target point above the 
Fischer e11ipso1d-ft 

IGUIDE-first pass indicator for E6RT 

XLFAC-Load factor in fps^ 

DELAZ-heading alignment error-deg 

HLS-height above landing site-ft 
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A-Array 

1197 

1198 

1199 

1200 
1201 
1202 
1203 

1204-1206 

1207 

1208 

1209 

1210 
1211 
1212 

1213 

1214 

1215 

1216 

1217 

1218 


Definition 

TRANGE- range to way point 2-NMI 

RDOTGN (ROOT) - h-dot with dispersions-fps 

2 

DRAGG-actual drag a.cce1eration-fps 

VTAEM - nominal velocity for TAEM initiation - fps 

MANP (AUTOPLAMP) - Manual pitch control indicator; 0-auto, 1 -manual 
MANRY (AUTORYLAMP) - Manual roll control indicator; li-auto, 1 -manual 
MANBF (RETRACTBF) - Body flap retract flag; O-active, 1 -retract 
Not used 

Velocity at v/hich to stop alpha-wind correlation-fps 
Not used 

NIC -commanded Nz change-g's 
Commanded angle of attack-deg 
Initial angle of attack-deg 
Navigated angle of attack-deg 
Not used 

K-gain on alpha rate 

2 

Ka-gain on alpha error (co N^) 
maximum alpha rate-dps 

niaX 

2 

alpha acceleration-dps 

fUaX 

2 

a-actual alpha acceleration-dps 
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Definition 


A-Array 

1219 

1220 
1221 
1222 
1223 
1224- 

1225 

1226 

1227 

1228 

1229 

1230 

1231 

1232 

1233 

1234 

1235 

1236 

1237 

1238 

1239 

1240 

1241 


a-actual alpha rate-dps 

Speedbrake deFlection in hinge line plane (0 to 98.6*^) deg 

Reference (nominal) SB def in hinge line plane-deg 

Commanded SB deflection in hinge line plane-deg 

SB deflection in waterline plane (0 to 87. 2°) -deg 

Landing gear flag; o-Up, 1-Down 

Body flap retract position-deg 

Body flap retract limit-deg 

Landing gear increment post-multiplier 

Landing gear deployment rate-1 /sec 

Not used 

Commanded bank angle-deg 
Initial bank angle-deg 

IWOW (WOWLON) - weight-on-wheels flag 
Not used 

K^-gain on bank angle rate 

2 

K^-gain on bank angle error 

A -maximum bank angle rate-dps 
^max 

2 

$ -maximum bank angle acceleration-dps 
^max 

2 

Actual bank angle acceleration-dps 
Actual bank angle rate-dps 
Not used 

Actual bank angle-rad 
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1248 

1249 

1250 

1251 


Definition 

DRPRED - range to minimum energy pt. - ft 
DNZCF ~ filtered unlimited - g’s 
EN - reference energy over weight - ft 
PHIC - unlimited roll command - deg. 
HERROR - alt. error - ft 
HREF - ref. altitude - ft 
QBERR - dynamic pressure error - psf 
QBREF “ ref dynamic pressure - psf 
Not used 

THETAN - navigated pitch angle - rad 
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A-Array 

Definition 

1252 

PHIBN-navi gated roll angle-ioad 

1253-1270 

Aero component table look-up results 

1253 

*k 

^•-BASIC ^f^BASIC 

1254 

%ASIC %ASIC 

1255 

'"'"BASIC 

1256 

AC|_ Speedbrake (SB) 

1257 

^'^Dsb 

1258 


1258 

ACmsB 

1259 

AC|_gp and ACj\|* Body Flap (BF) 

1260 

ACDgp and AC^gp 

1261 

ACrngF 

1262 

ACLpL 

1263 

A'Copj_ and ACA^p 

1264 

ACfOpL 

1265 

AC|_ Landing Gear (LG) 

1266 

^'^D-LG 

1267 

^Vlg 

1268 

•ACg "^Ground •’Effects (GE) 

1269 

^''D-GE 

1270 

^'^m-GE 


* Aero Data tape contains normal and axial force coefficients in the viscous 
intraction regime. TAERO does necessary conversions after look-up. 
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Definition 


A-Array 

1271 

1272 
1_2_73 

1274 

1275 

1276 

1277 

1278 

1279 

1280 
1281 

1282 

1283 

1284 

1285 

1286 

1287 

1288 

1289 ' 

1290 

1291 

1292 

1293 


X-cg position-inches aft of forward position 
Z-cg position-inches below nominal 
Scheduled eleven deflection-deg 

Elevon deflection used in aero components determination-deg 
Body flap deflection used in aero component determination-deg 
Height parameter for ground effects determination 

Mach number switch value for inviscid tables 

HREFTD (HREFTDOT) - ref. alt, rate-flare/shallow glide - fps 

HDTERR (HDOTERROR) - alt. rate error - fps 

HRFDOT (HREFDOT) - ref. alt. rate - final flare - fps 

Commaned bank angle-deg 

Flag for VECWIN routine: O-Off, 1-On 

Wind type selector 

Probability of vector wind ellipse 

Clocking angle from North-deg 

Shear/ Gust altitude-ft 

Probability of conditional shear circle 

Profile selector (geographic) 

Profile selector (month) 

Shear factor 
• Gust factor 

Delta clocking angle on shear circle 
Delta clocking angle on gust 
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A~Array 

1294 

1295 

1296-T299 

1300 

1301 

1302 

1303 

1304 

1305 

1306-1308 

1309 

1310 

1311 

1312 

1313-1317 

1318-1319 

1320-1322 

1323 

1324-1326 

1327 


Definition 

Magnitude of resultant wind velocity vector-fps 

Wind azimuth-deg from North in direction .wind is blowing. 

Not used 

3 

Previous value of atmospheric density-slugs/ft 

Previous value of altitude-ft 

d /d -s1ugs/ft^ 

P n 

dNj/dy-l/sec 
Gain on df^.j./dy 
Gain on Nj 

Interim calculations in open loop guidance 

Delta load factor-g's 

Actual angle of attack-deg 

Predicted load factor-g's 

Gain on predicted 1 o.ad -factor 

Interim calculations in open loop guidance 

Not u-sed 

Interim calculations in open loop guidance 
Gain for h-dot variation 
Interim calculations in open loop -guidance 
Gain for density variation 
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A-Array 
1328-1331 
1332-1339 
1 340-1 345 
1346-1348 

1349 

1350 

1351 

1352 

1353 

1354 

1355 

1356 

1357 


1358 

1359 

1360 

1361 

1362 

1363 

1364-1379 

1380 

1381-1389 

1390 

1391-1397 

1398-1399 


Definition 

Interim calculations for open loop guidance 
Not used 

Interim calculations for open loop guidance 
Not used 

2 

Commanded heat rate-BTO/ft sec 

Commanded load factor-g's 

Ambient temoerature - °K 
Temperature ratio 

Viscous interaction parameter - Va> 

Proportionality constant in temperature - viscosity relationship 
LAM 

L/D 

Not used 

Viscous interaction switch-over value for aero data 

Aero indicator for CLCT : O-CL^ CD from tables, 1-use TAERO 

Guidance mode indicator: -3-ADC, 4-TAEM, 5-A/L 

0-con » g, 1-con. (j;, 2- 4> from table 
Bank indicator in CLCT: 0-Bank right, 1-Bank left 

Cj^ total 

Load factor limit-g's 
Not used 

Actiral’ bank angTe-deg 
Not used 

Heat rate post-multiplier 
Not used 

Interim calculations for open loop guidance 
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A-Arra,y 

Definition 

1400 

Flag for TPS calculation; 0-no calculations, 1-do TPS 

1401-1433 

Surface temp,eratures-°F 

1434-1438 

Bare top temperatures-^'R 

1439 

Not used 

1440-1444 

Bare top temperatures-^R 

1445-1451 

Not used 

1452 

Transition Re-1 /ft 

1453 

Turbulent Re-ft 

1454 

Normal shock Re-T/ft 

1455 

Panel heat ratio 

1-456 

Turbulent heating multiplier- 

1457 

Not used 

1458-1461 

TPS stop conditions 

1462-1494 

•Maximum surface temperatures-°F 

1495-1499 

Maximum bare top temperatures-°R 

1500 

Not used 

1501-1505 

Maximum bare top temperatures-°R 

1506-1538 

Time of maximum temperature-sec 

.1539-1543 

Time of maximum bare top tempera ture-.sec 

1544 

Not used 
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A-Array 

Oefini tion 

1545-1549 

Time of maximum bare top temperature-sec 

1550-1574 

Panel LRSI weights-lbs 

1575-1599 

Panel HRSI weights-lbs 

1600 

Angle of attack-deg 

1601-1625. 

Panel thickness-in 

1626-1650 

Backface tempera,ture,s-'^F above 350 

1651 

Number of panels 

1652 

Panel area-ft^ 

1653 

Heat rate-BTU/ft^/sec 

1654 

Frees tream gas tempera ture-°F 

1655 

Panel number 

1656 

Panel thickness-in 

1657 

Normal shock viscosity 

1658 

Heat load-BTU/ft^ 

1659 

Total TPS weight-lbs 

1660 

Norma-1 shock pressure-atm 

1661 

Normal shock enthalpy 

1662 

Not used 

•1663 

Deflection factor mul<tipller 

1664 

Body flap deflect! on-deg 
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A-Array 

1665 

1666 
1657 
1668 

1669 

1670 

1671 

1672 

1673-1677 

1678 

1679-1683 

1684-1716 

1717-1749 

1750-1800 

1801-1802 

1803-1804 

1 805-1 806; - 

1807-1808 

1809-1810 

1811-1814 

181 5-1818 

1819-1822 

1823-1826 

1827-1828 

1829-T830 


Definition 

Eleven deflection-deg 

Deflection angle-deg 

Panel design heat load-BTU/ft^ 

Panel heat pulse time-sec 

2 

Total panel area-ft 
Average panel heat load-BTU/ft 
Total LRSI weight-lbs 
Total HRSI weight-lbs 
Bare top tempera tures-°F 
Not used 

Bare top temperatures-°F 

Panel heat loads-BTU/ft^ 

2 

Panel heat rate-BTU/ft -sec 
Not used 

CUBC3 (CUBICC3) - coefficients for HREF - 1/ft 
CUBC4 (CUBICC4) - coefficients for HREF - 1/ft^ 

DREMAX (DELREMAX) - delta range used to compute EMAX - ft. 
EDELNZ - delta E/W for energy limits - ft 
,EDRS - slope of E/W for S-turn 

EMEPCl - constant E/W used for MEP determination - ft 

EMEPC2 - slope of E/W used for MEP determination 

ENCl - constant E/W used for ref. energy - ft 

ENC2 - slope of E/W used for ref, energy 

EOWSPT - range for energy ref. index determination t ft 

ESI - constant E/W used for S-turn - ft 
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A-Array 

1831-1832 

1833-1834 

1835-1836 

1837-1838 

1839-1840 

1841 -.1842 

1843-1844 

1845-1846 

1847-1848 

1849-1850 

1851 

1852 

1853 

1854 

1855 

1856-1870 

1871-1872 

1873-1876 

1877-4878 

1879-1880 

1881-1884 

1885-1888 

1889-1892 

1893 

1894-1949 


Definition 

GAMSGS - Steep glide slope angle - deg 

HALI - Altitude of Minimum Energy Point - ft 

HFTC - Altitude of Nominal Entry Point - ft 

PBGC - Lower limit of initial ref. glide slope 

QBC2 - Slope of y'ef. dynamic pressure with range - p'sf/ft 

QBRLL - ref. dynamic pressure lovy^r limit - psf 

N 

QBRUL - ref. dynamic pressure upper limit - psf 

RMINST - Min. range to start S-turn - ft 

RNl - Range constant of energy computations - ft 

XA - Steep glide slope ground intercept - ft 

Time from start of phase for ground intercept change - sec 

Velocity for ground intercept change-- fps 

Altitude for ground intercept change - ft 

Range for ground intercept change - ft 

Ground intercept change flag 

Not used 

DSBCTD- Speedbrake deflection at touchdown - deg 

HK - Constant -G circle center altitude - ft 

SIGMA - Exponential distance - ft 

VREF - Ref. airspeed - fps 

XK - Constant -G circle center range - ft 

R - Constant -G circle radius - ft 

XEXP - Exponential capture "distance - ft. 

H6EAR - Gear deploym.ent alt. - ft 
Not used 
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1950 

1951 

1952 

1953 

1954 

1955 

1956-1998 

1999 


Flag for sonic boom output:' 0-no sonic boom; 
l-do‘ sonic boom 

•Mach # to start sonic boom output - phase 1 
Mach # to stop sonic boom output - phase 1 
Mach # to start sonic boom output - phase 2 
Mach # to stop sonic boom output. - phase 2 
A-array location to qua'nity 1>or output check 
Values of A(A(1955)) -at .which- .output -i-s to be done 
Counter for sonic boom output 
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DETERMINATION CONTROL SURFACE DEFLECTIONS 
AND TRIM AERODYNAMICS 


GET AERO FOR EL] AND BFn 


EL Eleven 
BF ■ Body pfap 


EL-j - VALUE OF ELEVON USED IN TABLES 
ELg - SCHEDULED ELEVON DEF. 

BF-, - VALUE OF THE BODY FLAP USED IN TABLES 


IF BF IS IN RETRACT MODE 


INCREMENT EL • 


IF BF WAS LIMITED 
LAST PASS 


INCREMENT BF 


GET AERO' FOR NEW BF-, 


TRIM CONDITION FOUNp' 


GET AERO FOR NEW. EL. 


IF TRIM CONDITION FOUND ^ > COMPUTE CLj, AND RETURN 


IF EL AT LIMIT 


GET CLOSEST AERO AND RETURN 


COMPUTE 6Ej, CLy, Coy, AND 


RETURN 


INCREMENT EL 


6 4. Trim eleven deflection 
et . 

(Sbft deflection 

Cd^ Trim drag coefficient 

Cl frim Tift c"^ficient 



















TABLES FOR ENTRY -TO-TO.UCHDOWN SIMULATION 


TABLE NO. 

QUANTITY 

INDEX 

1 

Wg(0)^.g, ft/sec 

335 

2 

W'^(0}jg, ft/sec 

336 

■5 

W^(Q)kg, ft/sec 

337 

4 

P / P' 

338 

5 

•P « /e<L 

339 

.6 

c /O' 

340 

7 

- Pco/ Pj 

341 

8 

P'^ , Ibs/in^ 

342 

9 

- slugs/ft^ 

344 , 

TO 

- ft/sec 

345 

n 

- lbs sec/ft^ 

349 

39 

K, Heating Option 1 

597 

40 

Scheduled elevon-deg 

1221 

41 

(j)^ - deg 

1380 

42 

- deg 

1210 

46- 

NTc - g's 

1363 

50 

Scheduled speedbrake -deg 

1273 

51 

<^Lbas1c - ■'™ 

1253 

52 

‘^Dbasfc - 

1254 

53 

'^mbasic - 

1255 

54 

C{ — Speed Brake (SB) - low mach 

1256 

55 

Cn - “low mach 
SB 

1257 

56 

CmsB " mach 

1258 

57 

— Body Flap (BF) - low mach 

1259 

58 

Cpgp - low mach 

^^®°ENCL0SURE 4 
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TABLE NO, 

QUANTITY 

INDEX 

59 

Cmgp - low mach 

1261 

60 

Cl — Elevon (EL) ~ low mach 

1262 

61 

Co mach 

1263 

62 

Cju PL" mach 

1264 

63 

Cl Landing Gear (LG) 

1265 

64 


1266 

65 

LG 

1267 

66 

Cl — Ground Effects (GE) 

1268 

67 

‘•'’ge 

1269 

68 

CrtiQE 

1270 

69 

'^'-basic ■ 

1253 

70 

'^“baslc - •''9'' “9'’ 

1254 

71 

Hasic - '’’9'' ”99'’ 

1255 

72 

CLjg - high mach 

1256 

73 

Cg^g - high mach 

1257 

74 

Cm^g - high mach 

1258 

75 

CLgp - high mach 

1259 

76 

Cogp - high mach 

1260 

77 

Cmgp " high mach 

1261 

78 

CLp^ - high mach 

1262 

79 

Cqel ~ high mach 

1263 

80 

CmpL " high mach 

1264 

81 

Cq inviscid uncertainty 

1061 

82 

Cl inviscid uncertainty 

1062 
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TABLE NO. 

QUANTITY 

INDEX 

83 

Cd viscous uncertainty 

1061 

84 

Cl viscous uncertainty 

1062 

85 

altitude rate error-fps 

1055 

86 

Range-to-go - NMI 

1057 

87 

altitude error - ft 

1041 

91 

'"basic 

1253 

92 

%asic 

1254 

93 

'"basic 

T255 

94 

C|\|gP viscous 

1259 

95 

CAgp viscous 

1260 

96 

Cnigp viscous 

1-261 

97 

C^ipj^ viscous 

1262 

98 

Capl viscous 

1263 

99 

CniEL viscous 

1264 

101 

Panel design heat loads 

1667 

102 

Normal shock viscosity 

1657 

103 

Panel heat ratios 

1455 

104 

Panel Reynolds no. 

1452 

105 

Panel thickness - in 

1656 

106 

2 

Panel areas - ft 

1652 

120 

CRsb - low mach 

1094 

121 

CHgp “ *^3ch 

1091 

122 

Cgpj - low mach 

1092 

123 

CHpo " 

1093 
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TABLE NQ. 

QUANTITY 

INDEX 

124 

“ high mach 

1094 

125 

Cfigp - high mach 

1091 

126 

Ch^j - high mach 

1092 

127 

CHpQ " high mach 

1093 . 

The following 

are required if CbCT ts-used >and-A'( 1-3595 

= 0 

16 

C Q - i n V i s c td d rag -co ef f i .ci:e rrt 

■417 

20 

C|_ - inviscid lift coefficient 

1362 

44 

Body Flap Schedule 

1664 

45 

6^^ " T^iviscid eleven deflection 

1 665 

47 

C[) - viscous drag coefficient 

. 417 

48 

C|_ “ viscous lift coefficient 

1362 

49 

6^1^ - viscous eleven deflection 

1665 
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Z RUN M70663,EXA3,BETA,3304J,E062 ,P,10,.6 



A5G A=X07475 
ASG F=XiA728 
— ASG H=X17990” 
- ASG C=C 

ASG-tS,!3tG,X,Z- 

XGT-CTjR 

TRW A 
IN .’ A 



XGT SGGEM 


X~1'4 0 0 0T2 0 OlTl' S^TIT' 


#|- 


F 133 1000, 1231000, 200012 U 2322213 ,12223134 
FI 132000, 1323212,2223312 a 1330004 

— F 200 0231, -2 000221, -131 1000,-200012 14 : ^ 

F1131000,0, 2000331, 20-00321, 2C002114- 
F2000231, 2000121,2000331,20002214 ' '' 

F1131000, 2000321,1311000,0,20002114 ■ 

F'0r3'7141 59 27723 0072 Q 0'071T57 079 £324“ 

F 133 1000, 1.321000,2000211,20001214 
F 1121000,0,1231000, 12210 00, 131 1 0 004 
F X.31I000, 1231000,1221000,200021 14 

-^F i33 100 07132LOOOVO, 2 000121,-ri310 004"-~' ~ 

F 13 3 100 0, 13 210 00,20 00 211,2232 12 3 ,11332 22* 
F 131 20 00, 1122 233,2 223.132, 13130004 
FI 133000, 1223 3 12, 2212323, 2000 132 ,23132224 

~-nn 2 2zmn z ro o 072 0 0 o 237 ti 3 - 3 to 0:0-4 — ^ 

F 1,4, 1311000,0, 2000211,04 
Fill 1000,0,1211000, 0,l'3ll0004 
FI200,20,38,25,164 

~'F 4 ;-0 9 2 2 3 ,"0757023 94"rO-6T0737“4-4-03-2-9E'-=34-T — - 
F 0 ',c,' 0 , 0 , 0 , 0 , 0 , 0 , 0 , 0 , oi 
. F2Q925741, 20855591 ,7, 2921-15E-54 : ',7 

F32,146 502, .05218839,4. r60188E-54 . - •' 

— F'2 09 029T07327ZT."7 75’07 67T033757729"57794--- 

F 14 ,695833, 2.37692E-3, L 1 16. 4 , 3 . 737 17E- 74 
5 F32.174S 


X47l$ . 

X1,5$ 

"FI"*, ar,3-,"i* 

FO's 4000,5000,8000, 


ru, iuuu,;,uuu,ouuu, xu 

F65, 65, ,98. 6, 98. 6, 03 

"x-4'a$‘“" 

X17J E 

F 1,0,6, 14 

FO.g, 1.5,3,5,8,12,04 
■•F4r7"^7-37275*,‘2".'5T=^2.'ir 


5PEEDSRAKE SCFEDUIE, 

10003 ^Q4 


ELEVON SCHEDULE 



CT22-26 

CT27-30 

CT31-34'-" : — 

CT35-39 • , ' ' ' 

CT40-43 . 

CT44-48 '7 • ■ ■ . 

"C T49-53— ro ,-PT7NrA-BTATP r/'2T~ 
CT54-57 . 

CT58-62 

CT63-66 

-C T67 -7 1"~— 

CT72-76 . ■ • • 

CT7 7-80 , . ■ 

CT81-85 

■CIBG-SO’' : 

CT90-94 
CT 95-99 

CTlOO-104 (X,Y PLOTS) 

CT105-1 09 — 

CTliO-120 (NOT USED) 

CTL-3 (A,B FOR FE, OMEGA E. 
CT4-6. { A0,A1', A2 FOR^FE)', 
T,T7"-10 ( RE'i G7FT/‘NMi'DEG/RA0 )“ 
CTl 1-1-4 (PSLjRHOSLtCLS-^’MUSL) 
CT15iGWM,CTl,6-^2L COM.-INT.*) 


PLOTS) 
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Xi ,0,0r 1$ . ■ * , ' . 

H NT R'V^ il-nnu CKDOW M~r-tMUrATT D N'T'A D C~G U I D A N CT' 
X0'vi<32t55, U3, 1,33$ 

"F'ir8i'55? — — : — — 

R.50,i’!= ■ ; . • - • 

F5I, 0« 

F52, 30* • . . ' 

T53 ri';9'2* 

F56,'0* 

F5 7,,0* 

F5S, 0* 

"F 59 , 3 : : — ”■ 

F6 3, 1* ^ 

F64,2* ■ ' , 

F65,7» , 

FLOOtO, 101 ,0* 

F102,0* 

F106,rlT6.83, 107, 2 6 . 1 5 , US', 40045 1* 

T 1267 ?. 5 549 ’,'IB 0 ',~ 1 T '4 8 ? raiT- 62'75 3 *"-- "T"— : 

F 12 U, 40, 1231^0* . ; 

F 1664,5* . .. . • 


Fli6a,2* 
F,1358,,005* 
F1277, 1«2* 
T12 2 67^*3 
F 1228,. 1* 


F 6 - 5.0 
! :~F70X 

-F 7 G 7 

F 7 10 

B|~TTn 

Pill 

[@■'=0 — F728 
8 ^ F 73 1 
PG F 735 

p 

F 745 

F 749 

- P752 

I F756 

I F 759 

' F 76 3 


, 1,653,200,656,2* 

76 ,702 r03'^< iTT. 

,0100,704, 0106, 705,* Q120, 
,0127:>,70 8,13i0,709, 1,210* 
, 1170,711,0110,712 7 0121 , 
70 nm r5TT3'8 0T“71 6'TT2 8 1 *■ 
.0106,718-, 0107,719 7 0122, 
,0128, 722, 1665,723, 1664* 
,0401,725,0109,726,0418, 
, 0 6 1'4T7 2 9 iTC 8 2 77 3 CT i;C 8 1 * 
,0402,732, 1353,7337 04217 
, 13 11, -736, 1085,73 7,- 122 0* 
,0591,739,0411,740,1343, 
rn'79T7“43 rro 9*6T74'4Tr22T*' 
,05 92, 7 4‘6^ 142 6, 7 47,' 1431 , 
, 1184, 7 50U 199 ,751,* 1313* 
, 14 29,753, 1428,754,1432, 
,"1 1"8 1'7~7 5 7 7 1 1 8 2 77 5 8 71 1 8 3 * 
, 142 7, 760,1430,76 1, 1433 , 
,109 9, 7 64, 0 579,76 5,1197* 


706 , 0117 * 
713 , 018 0 * 

720 , ciei* 

727, 1362* 

134,0417* 


748 ', 

755, 


1356* 
0351* 
13_5_L^ 
10 9 8* 


r*’ A“D ir."^ E Q UAT^rN T-E'G R AT EO 

OBLATE EARTH ' 

WIND/ATMO INDICATOR 

' 62 US STANDARD ATiMO. 


RUNGE-KUTTA INTEGRATION 
CONSTANT TIME STEP, 

SI NGLE PREC ISION 

— highest- OUTPUT' ROUTINE"” 

AERO LOADS COMPUTED ' ' 

. / ■ PUS. INPUT UPT ION,' GEOD. 

. , ■ VEL INP UT-'OP T IN ER 4 GD , 

CURRENT £ I N I TI AL”flMES 
RELATIVE LONGITUDE, 
INITIAL .LONG., LATi', ALT, 

INITIAL" VELOCITY VECIOR' 

■' ■ ; INITIAL ALPHA AND:;.PH 1 

INITIAL BF-PUSITI'ON 

I CLr- WET'GHT 

MID. C.G. 

REF AREA 8 LENGTH . 

POSTfMULTj PL liRS , 

QB'AR TO START TRI^ A>‘ERO 
VBAR SWITCH ■ 

MACH FOR AFRO TABLE SW. 

P j 

LG DEPLOYMENT RAT;E 
. OUTl',OUT 2 ,QUT 3 fMeQ. 

"-nDUTf^U'T — ■' r— 

TifeVAL'T.ViREL) ,V(A 1 R) ■ , • 

V,UN) , ALPHA! ACT)', ALPHA! CO «4 
■' FH 4 SE,R( EC) ,GAM(RGD) ,GAM(;iGC) 

— “G A 7 W I GD JVPH n ACT ') '7 PHT ! CGrrff 

LO(^G ,LAT!GD 1 ,AZtRGD) ,AZ! IGC) 
AZ( 16 D),DELE,DELBF' 

M,*LAT!GC) ,CY, CL - 

— NTiHMIE IHMBP ' T" 

• . QDAR rVBAR,e Z,C 0 !CX) •• 

■ Ni;!E ST] ,HMOE,DELSB!HL) 

, / QO.ur ,REtRHOV!M) ,LID 

LDOV ,HMSB,DELS.B! WL 7 = — 

QTUCPl ,TCP6, RHO 
ALDR EF ,ORAC iHDDT 
TCP4 ,TCP3,TCP7,TAM3 
T~‘D23,DREF,HDREF ”■ — -r— • -r~ 

' ! TCP2 ,TCP5 , TCP6,DP ■ r 
CPV, R!GC) ,R!TG) ' ’ 
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V 


F1236,5i 1237, 1.7=5= , 

FL216, 5, 1217,1* 

n= i 2 1 47 2‘ ; 1'2 IT, "2* 

F 12 3 4,. 7365, 1235,. .27 12* 


F 1360, : 
T17 3,“-: 
FlLOl', ■ 
FU027' 
FI L03,- 

“Friar, 

F'l L 4 1 , : 
■F-U42.,- 
F115'6,: 
“F.U47,; 
F1149,. 
F1150,l 
F a 5 2 , , 
"FclTT 47 ■ 
F'l I 57, , 
F1174,l 
F1175,; 
'F 1176,-: 
FIL77,: 
F1L69,( 
F11707- 
-FTI71T; 
F1172,; 
F1L92,; 


117.'8‘4T17'4' 
4.6822224, 
T*??. 337005 
40,1113,0, 
10200Trr3'2 
3 * 

-2'.4ai88, I 
5 000, ll^f4, 
20000, 1148 
31* 

0 * 

21, 1153,20 
9 500nT5T7 
.02, 1158, . 
80* 

2 * 

£ ^ — „ 

16.5* 

0 * 

-169.84, 11 

'20000* 

31993* 

2220 * 


1111,3.153 3738-3 m 21, 0* •' 
i, 11 12, 0.0192356, 1122, -7.88341 
1123,0* 

17 1’22 0 o“,-rr33727oo 0**“; ; — " 

143,-2',99933* . V , , 

3053 8.5,114 5,2383,1 146,30538’, 
1^2588 573 V 1T7 3 ,23 C 0 C n2 0 0 , 2 50 l 

1.96,115^2.086* • 

■57“0’602C7* : 

02, 1159 a 0 25 a 160,. 03* 


80,17* 


~T 1“40 07T*'"“ : — — 

F59 0, 2* ■ =■ ■ ■ . 

F 59 5, 1.0* ' ■ , , ' 

F1390,.82 . ■ ■ ■ • 

Tr4 3-4T<7a OTT435TTrOTT42r(rr6irDTr4 JTT6 Wfl'43 0 76-Er 0 *“ 
F 14 4 0,680, 1441,680,1442,680,1443 ,680,1444,68 0* 
F1458,591,- 1459,108,1460,-1 .0 ,1461,-8000 0* 

F 1 6 5 1 ♦ 2 5* ' J 

F 10 8 0, I* - . ■ ' 

F 10'88, 10.935, 1089, 19047, 1095,9331* BF,EL,SB ARE/ 
F1Q90,4* • , ^ 

F 13 12,1.3* 

F1323,-.0l00, 132 7,2. E6, 1304,-1250,1305 ,.5* 
8600, 1,601,100,602,4$ _____ 



MAX ROLL RATE & ACCEL 
MAX ALPHA RATE 8 ACG6L 

'•■'KAL7HAD0T7KALPHAERR0R 

KPtiIDOT,Kf^HIERROR- 

GUIDANCE Mode ind-icator 

- TARGET -LONG -a’^'L AT ^ — 

CALP0,CALP1,CALP2~QUAD1 
CALPC1,CALP1,CALP2-QUAD2- 
CALPQ,CALP1 ,CALP2-QUAD3 

“VA'LP"( 1-3 ]■"""" 

NALP • 

AKVAKl 

VQ'i VA',VAl,VA2 .i. 

VB1,VS1',VRDT,VTA-EM ^ 

ALFM 

PR'EBNK 

OF ’^RPTl ,EEF4 '' 

~ VTRANvETRAN =" — 

GS1,GS2,GS3 ,GS4 

ACNl . • ' 

DO L I M 

02 30 '• 

XWP2 ■ . , 

' RW 17 AZIMUTH AND ID 

“RTURN ^ 7='’ = 

DBAR 

EAFB HT ABOVE FE 

“"'TP'S' CAL'CULTiTiaNS ^ T”'’ 

HEATING INOICATOR 
VEHICLE NOSE RADIUS 
HEAT: RATE MULTIPLIER ' 

”1 N TTT AL” BAR E“T 0 P = 

TEMPERATURES 

TPS STOP CONDITIONS 

^ :.P ane ls _•_ 

CALCULATE HINGE MOMENTS. 

-ARMS FT2-IN- . ' ■ ' 

D ELTA EL FDR AILERONS ~ , ~ ■ 

GPRED. AMP. FACTOR 
QOOT 8 G CONTROL GAINS 
PHASE STOP CONDITION 

r~rr~T o 

. dj 

• . ■ , . > V ■ ' ' ; to 

- ' rp 

I > ■. ^ CO 

o 

“h 


CO 
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l4U?r (Utv*M<si r nrmt, *nc 


J< JLjl. 2 ^ 01 

1182 REF CRAG-FT/SEC 2 
U 99 ACT CRAG-FT/SEC 2 

1183 PEF HDGT - FPS 

•119 8 “”ACT~KD 0 T’~-"FPS 

1179 LODV 

1181 023 

1'19 7 RANGE-TO-GO-NMI 

" 13 ' 56 '“L 7 D 

1210 ' ALPHA COM - DEG 
;;100 TIME - SEC 
120 ■ V RELAT IVE FT/SEC 
" ro a — A L T I T U D E “ -“'F T ' 

'402 DYN PRESS LB/PT 2 ' 
166 ^ RF DEFLECTION- DEG 
1665 EL DEFLECTION- DEG 

131 0 ALT H A^—DEG 

128 1 • RANK ANGLE - DEG 
1380 RANK ANGLE - DEG 
614 LOAD FACTOR - G 

- 1 2 1 — G AMM A“=" D E G 

L 081 HINGF MOM BF IN-rLB 
1082 - IN EL HM - IN-LBS 
1083 OUT EL HM - IN-LDS 
~ 5'9 X — Q-DO r-BTir/ FT -2 / S EC— 


r0691T~l“ 

r0i08,~2 

,0120?+! 

» 0402^-1 
ro 1217=1“ 
,01G8r-l 
t 1665,-2 
, 16 64, + 1 
a 1 977-1” 
, 1356,-1 
,1179,0 
,1182,-1- 

71X9970 

, 1183,-1 
P 1198,0 
, U 81,-1 
7 V 1 1027 - 1 “ 
, 1081,-1 
, 1082,-1 
1083,-1 
r 13107^.“ 
, 12-10, +0 
, 1380,-1 
,12Bl,+0 
, 06147 ’- 1 $“ 
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1482T ' [lmiii«‘w Fwfcu Toe. ^ 


F40,04rOl,0120;00G0t 
^50, O'4t'02,O‘4OT,OO00V 
F5i, 10,36,0401, 1310, 
F52, 10,37,0401,1310, 
F53, 10, 38,0401, 1310, 

F54 1071 223', 

F55, 10, 40,1310, 1223, 
F56, 10,41, 1210, 1223', 
F57, 10,42,1310, 1275, 
F5e/i0,'43;i3l0,‘1275,' 
F59, 10,44, 1310, 1275, 
F60,10, 45,1310, 1274, 
F61, 10, 46, 1310, 1274', 
F62, 10, 47', 1310,1274, 
F63 , 10, S2, 0410, 1310, 
F64,.10, 53,0410,1310, 
F65, 10,54,0410, 1310, 
F6'6,10,~55', 131071276', 
F67, iO, 56,1310,1276, 
F68,.10, 57, 1310, 1276, 
F69, 10, 24,0401, 1310, 

F^O, 10, 25,-04Q1-,-1310,' 

F71, 10, 26,040 1, 1310, 
F72, 10, 27, 1310, 1223, 
F 73, 10, 28,1310,1223, 

F'74",710 ,“29T13T0TT 2 23T 

F75, 10, 30, 1310, 1275, 
F76 , 10, 3i,.1310, 1275, 
F77,10, 32, 1210, 1275, 

FT8 ,1 0, -33713 1 0712 7 Tr 

F79,10, 34,1310,1274, 
F30, 10, 35,1310, 127/*, 

■ F$ )l, 10,01,1310, L353 , 
F9-2'n070‘2Tr3T07l35'37 
F93, 10,03,1310, 1353, 
F94, 10, 04,1275,1353, 
F95 ,10,05,1275, 1353, 

F96 , 1 0 ,-06 ,' l 275713-53T 

F97, 10, 07, 1274, 1353, 
F9S, 10, 08, 1274, 1353 , 
F99 ,10,09,1274, 1353, 
Fitnrro qitittb55to o oo’ 
F 102,08,18,1660-, 1661 
F103, 03, 19, 1655, 1600 
F 104,00,20, 1655,1600 

F10 5,03,27,"I'6'55“,r6 58 

F106, 00, 28, 1655 ,0000 
F120, 10,48,1310,1223 
F12i, 10,49, 1310,1664 

FITI TOTS 0 1"6 (T5" 

F 12 3, 10, 5 1,13 10, 166 5 
F 124, 10,20, 1310,1223 
F125, 10,21, 1310,1664 

F126, 10', 22, 1310,1665 

F127, 10,23, 1310,1665 


0000,1221 ,00 ,058 7,00,000 0 
"00007127 3", 0 0 ','0 5 87 70 0 T 0 00 0 
0000,1253,00 ,0587,00,040 1 
0 000, -12 54, 00 ,05 8 7, 0 0,0401 
>0000,1255,00,0587, 00,040 1 
. 0 4 0 171 2'5 6 ,'0 O' ,'0 5 87 70 0 , 040 1 
>0401,1257,00 ,0587,00,0401 
0401,1258,00,0587,00,0401 
C40r, 1259, 00, 0587, 00, 040 1, 
040 171 2 6 C70 0 ;■ G 5 87", 0 0 , 040 1 
04 01, 1 2 61, 00, 05 8 7, 00 j 040 1 
0401 ,1262*00 ,0587, 00,0401 
>0401,1263,00 ,‘0 5 87,0 0,0401 
0401, 1264 iOO ; 0 587, "00, 0401. 
0000,1265 ,00 ,122 7,0 0,000 0 
COOO, 126 6, CO ,1227,0 0,0 00 0 


12 i4, 1268,00 ,0587r00,,0000 
127 4‘, 1269,00 ,05 87,0 0,000-0 
127 4,12 70 ,00 ,05 87,00,'0000 
CCOO, 12 53,0 0,0 587,0 0, 1363 


0401 f 12 56 ,00,05 8 7,00, 000 0 
0401712 5 7,00 ,05 8 7,-00,000 0 
‘040 1‘7X'2 5870 OTO 5 8 770 0 ,'0 00 0 
G4Ql;i259,00,0587,00,+35 3.- 
040''1‘;1260,00 ,05 3 0 0,-135 3 

0401,1261,00,0587,00, 1353 
0 40 1-7 12 6 27 007 C 5 87 i'O 0 , 135 3 


040 1.- 12 63 , 00 , 0 5 8 7, 0 0 135 3 
C40 IS 1264,00 ,0587,0 0*1353 
CCOOS 12 53,0 0,05 87, 00; OO'O'O 
G'CO 0TOT470'D"70‘5 8 77O"0 ,‘p0'0‘ 
000011255,00 ,05 87, 00, '0000 
131; 0412 59, 00, 0 5 87,0*0*000 0 
13105 1260,00,05 8 7,00, 000 0 


:oo'oo 
Oooo 
00 00 
00 00 


0000 , 000 $ 
0000 , 000 $' 
1277, 069$ 
1277,070$ 
1277,071$ 
1277,072$^ 
.1277,07 3.$ 
>1277,074$ • 
1277, 075$ 
1277,076$" 
T277, 077$ 
1277,078$ 
1277, 079$ 
1277, 080$ ' 
0000 , 000 $ 
0000 , 000 $ 
:> 0000 , 000 $ 
0000 , 000 $ ■ 
0000 , 000 $ 
'0000, 000"$ 
1358,091 $ 
13-50, 092$" 
1358,093$ 
0000 , 000 $ 
0000 , 000 $ 

0 000 , 000 ,$“' 
1358,094$ 
1358, 095$ 
13 58,096$ ' 
1358-, 097$ “ 
135 8, '098$ 
13587099$ 
0000 , 000 $ 

0 000 ,“000 'I""' 

0000,000$ i 

0000, poo$ 
0000, 000$ 
000070001$ ‘I 
0000,000,$ 
0000, '000$ 

0000. 000$ I 
i7OOG0V0Q0$" 
,0000 ,000$. 
',0000,000$ 
,0000 , 000 $ 
, 0000 , 000 $ 

1, 0000 ,000$ 
,1277 ,124$ 

, 1277, 125$ 
,*12777126$' 
, 1277,127.$ 
',0000 *000$ 
, 0000 , 000 .$ 

0000 ,00.0$ 
, 0000 , 000 $ 


SB SCHEDULE-- 
EL"'‘SCHEDUl.E"'~ 
CLB LO MACH 
COB LO MACH 
CMB LO MACH 
■CLSB'"'L'0‘*"MAC'H*" 
CDSB LO MACH 
CMSB LO' MACH 
CLBF LO MAGH 
“CDBF‘~LO ‘MACH"* 
CMBF LO MACH 
CL-EL LO MACH 
CDEL -LO MACH 
CMEr"L0‘ MAGH'* 
CLLG LO MACH ■ 
COLG LO MACH 
CMLG LO MACH 
CLGET“LO “MACH” 
CDGE LO MA.GH 
CMGE LO MAGH 
CLB HI MACH 
COB — H r M'ACH' 
CMB HI MAGH 
CL:SB .HI, MAGH 
CDSB HI MAtH 
CMBF~Hr*MAGH~ 

CLBF HI m,ax:h 

COBF HI M'AGH 
CMBF HI MACH 
cL'eL-"Hr-MAr;H** 
CDEL HI MACH 
CM EL HI MACH 
CNB VISC-CIUS 
CAB-^YrSCEiSS"' 
CMB VISCDUS 
CNBF VISCOUS 
CABF VISCQOS 
CMRF-^VISGOl|S"~ 
CNEL VlS-COiS 
GAEL .VISCOUS 
CMEL VISCOUS 


■ NOR. . SHOCK VIS, 
HEAT RATIOS 
-TRANS. RE NO. 
THICKNESS"' — 
PANEL AREAS 
CHS 8 LO MACH 
'CHBF to MACH . 

"CHEr-Lo- mew . — 

CHEO LO MACH 
CHSB HI M'ACH- 
CHBF HI M;ACH 
"CHEr HI MACH “ 
CHFO hi MACH 
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■^TTTt'O'TT'S 7 

. r ENTRY-TD-TOUChDDWN SI I^ULATIQN , 7/lEM GUIDAMCE 
XQs, I2a, OrOfOf 33$ 


"FI 80 1 ; -8 rE^7T8 O'ZT'O’’/*" 
Fia03,4E-13Ti804,0’!' 
71 80 5»5A00Ctia06,O=? 

F 1 80 174000, 1800, 0^^ 

-c~i an c . "ACiQ/x A-1 av RT-n~ 


■TI823', r6005,.1824;0,l( 
Fifi27, 110000,1828,0* 
81828,90000, 1330,0* • 
Fig 3 1,-24, 1832-, 0* 

~F“1 8 3 3', 1 0 0 1 8 n 83 4 TO 

t: ^ a a lonin vo'^c. ni'- 


0 * 

1818,0* 

0 * 


-“~cmc3- 
eUBC4 
DREMAX 
• EDELNZ 
—EORS"’- 
EMEPCl 
,EMEPC2 
‘ENCL 
■”‘ENC2 
EOWSPT 
ESI ■ 
GAMSGS 
-■'HA'U I,“T 
WFTC . 
PBGC ■ 

, -,QP;C2 
":OBRLC~' 
■ QBRUL > 
RAUNST 
RNJl’ 

-^A- 




( 1 , 2 ) 

( 1 , 2 ) , • . - 

(1,2) ■ ■ 

“(1,21 — ^ = 

(.1,1) (2,1) ( i,2)-(2,2J 
(1,1) (2,1) Uy2J (2,2.) 

(1.1) (2,1) (I, '2 J (2,2 ) 

•( l,l) '(2,-l) Cr ,2J r2,2J- 

U,2) 

(1.2) 

( 1 , 2 ) 

“(142)“““ 

( 1‘, 2 ) 

( 1 , 2 ) 

( 1 , 2 ) 

'(1,2) r- 

( 1 , 2 ) . 

( 1 , 2 ) 

(l-,2) . . 

“( l, 2 )-r — rr 


r iv-ciyt — c o , i w w , u''- , i ^ f 

! 'F.ia«1 1,210) 1842V 0* ^ ^ ^ ^ j_^ 2 ) r: 

Fla^2, 300, 1844,0* • ‘ ’ ',0BRUL> (1,2) . 

FISTS', 152000, 1846,0* .RAtlNST (1,2) 

Fia<7,3(;456,6, 1840,0*. RNl’ (l-,2) 

“F-18 4gV~5000VT8SOTO* "Ti^A- “(1,2)— . t 

81851,2000, 1852,0, 1853, GU854,-2E6* . GLIDE SLOPE CHANGE CO^iO; 

770.1,6,702,63* ’,.1 ■ ' -.foUTPUT 

. “-F?03,aiGOV704','Or08?7GrrOl2’Or706r01T7*“““TTHE , ALT, V OREL) ?V(AI R)““iry7 

i. F7C7,.0126, 70S, 1310, 709.(rl2I 0“^ ' V.( IN ) i A L PHA ( ACT ) , AL PHA ( C Qf;l3 

' ' FTiO, 1178, 711, 0110, 71'2r0l21, 713, 0180* ' PIhA'S PiR ( EC),,GAM(,RG0 J ,GAMU'GC) 

^ . E714VG127,7I5rl380v7I'6:V128i* G I GO ) , PH I ( ACT 3 , PH I( C08)::. 

I F717 /QrQ6T7r3r0107;'7i9F01'22r7207TlUlF LTQNG V LAT (GD,)-i:AZ( RGD )TA'Z CIGC-)~ 

F 72 1,0120, 722vl'665, 723, 1664* AZ7 1 GD) , DEL E , DELB F ■ 

F72^, 06 01, 725, 0.109, 7 26, GM8,?27 , 13,6 2* Ha’AT(GC) ,C.Y,CL 

F 728, 0614, 729, 1082, 73CUC81* ' NI',H MIE , HMB.F' 

1 ,-0402V732,1244,-733V0^2-!:,''?34T06T7* QBAR , EN V CZV.'CD FCX ) “ 

F735, 1311,736, 1085,731, 1220* NI'(F ST) , HMOE , DELS 8 ( HL ) 

! F738, 1249, 739,0411)760) 1343,741,1356* ■ QblRE FiRETRHOV.’.lN), LID 

I F742, 1222, 743, 1096, 744VI223* . . DS;BG AT , HNS8;, DELS 0 ( WU 

^F74 5 71 2 48 ,'7 46-71-247 174771 246 7T48T03 51'* eB.ERR V H REF , H E R RO R7RH0 

• F749, 1209, 750, 1199,751, B13* IMZ'C, DRAG ,HD;0T . 

F 752, 1186, 753, 1187,7 54 , 1 18 8,7 55, 135 1* XRW, YR W ,ZRW ,7 A MB 
F756, 1243) 757, 1245,758,1242* . QN7.CF ,PHIC, DRPRED 

1 F-7 59 ,“U 89'776C7T190 ',76"1 7119 IT7624T0 FS**— ~XP:OTRWVYPOrR W , ZDOTRWTOP 

I F763, 1099, 764, 0579, 765 ,1197* CPV> R (GC).,R ( TG) - 

I. F1234, 1-1312,1235, .64* ' KPHI.DDT , KPHl ER 


C.Y , CL 

R P' 


F 1234, 1-1312,1235, .64* 

'nZ‘1 6V5rr21 772“0 *~ 

F1236,20,1237,50* 

•F53,-96* 

F590,0* • . ' 

■"F 6 024 5$ 


KPHI.DQTrKPHlERROR , 

A L”.PH A R A T E'~£ ~A'C G EC A X" “ 
ROLL RATE L ACCEL MAX 
INTEGRATION TIME, STEP 
STOP QDOT CALULATJON 
phase 'ST0P''C0NDITI0N ■ — 
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xr, 34»0$ 


U97 RANGE-TO-GO-NMI 
1210 ALPHA CGM - CEG , 

100 TIME - SEC 
-120-" V REL"ATTVE— FTTSEC" 
iOS ALTITUDE - FT 
402 ■ OYN PRESS LB/FT2 . 
1664 DF DEFLECTION- DEG 
“1'665“^- EL”DEFLECTTON-~'DFG~ 
1310 ALPHA - DEG 
120 1- BANK ANGLE - DEG 
1380/ SANK ANGLE - DEG 

-1209“ NOR-.LOAD’ COtM~G 

,614 ■LOAD FACTOR - G 
121 GAMMA - DEG 

108 1: hinge mom bf in-ls 

TO 8 2 1 N — EL”HM-'~TN~1TB‘S- 

, 1083 OUT EL HM - IN-LBS 
1187’ CPOSSRANGE - FT 
1 186 ■ DOWNRANGE - FT 

T183“H ABOVE RW~FT-““ 

1247 HREF - FT 
1223 SB OEFLECTION-DEG 
109 6, SB HM 4- IN-LBS 
T24 9 Q E A-R-REP'— P S"F~“ 


X‘9v 2 5S 


O O 


“X-,10 07 
5C10-0, 
X.12 0,, 
Xi20r 
"XIOQT 
X120, 
X120, 
Xl20t 
-X'1-20 ■; 
X 120 , 

XI 2 O 7 

X 1-20 , 
Ta'ZOT 
X12Q, 
XLZOr 
X 120, 
"X'120f 
X 120 , 
X120, 
XlOO, 
“’/TIOCV 
K 120 y 
X 120 , 
X1136 
“X"120, 


0 r08 7-2“ 
0120, +1 
402,-1 

1249. 0 • 
'or2r,“-i“- 
01G8r-.l 
1665,-2 
1664, -M 
1C81T-1"~ 
1082,-1 
1003,-1 
1310,-1 
■12T0T+0" 
1380,-1 
1281, +0 
1209,-1 
ll88,-r~ 

1247.0 
1197,-1 
U 87,-1 
ir86","0"^ 
1223,-1 
1096,-1 
,1187,-1 
06 14T-T$‘ 
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XITITOTr^^ 

entry-td-touchdgwn 

X0;C93,0.0,0.0,33$ 


Slf^ULATIGNt AUTOLAND GUIDANCE 


"Fnr7'T7 0 71 e 727“0*' 

F187 3, L7000a8?47 0t 1£75,0, 18 76, 

F 1377, 863. A, 1878, 0« , 

Cl n ", Q . c: n -3 . 1 u 13 n . rati. 


1 A.WW.'T ^ ^ 

P 1893,900* 


■"DSBCTD" 

:HK 

SIGMA 

VREF 

"'XK ^ 

R 

•XEXP 

HGEAR 


'(TV2') — : — - 

a-,1) (2,U (1,2) (2,2) 

11»2) 

( 1 , 2 ) 

- ( 1 , n (2Tir(l'T2-)7 272“)- 

(1, U (2,1M1,2) (2,2) 

(la) (2, 1) (1,2) (2i2) 


F 701 
F703 
F7Q7 
— F-710 
F714 
F717 
72 1 
"~F‘7 '2 
P7.Z3 
F7 3 1 
F7:3 5 

F73S 

F 7/t 2 
F7.45 
F7/,9 
_ 2' 
F756 
F739 
■ F763 


01 ,;6,702,63* 
03,0100, 7C4, 
07,0126, 7G8, 
10-7ir78",-7TlT 
14,0127,715, 
17,0106,718, 
21,0128,722, 
2<i-7C40If7-2'5,' 
23,0614,72.9, 
3 1 ,0402, 732, 
35, 1311, 736, 
3 3 -,-127 97739 7 
42- 12 22 f. 743, 
45, 1248, -746 ,• 
49,1209,750, 
5 27 U" 8 6", "753 7 
56, 1224,757, 
39 , 1189,760, 
63,1099, 764, 


01C8,7C5'i 
1310,709;,' 
•017-077 127 
1 380,716,'-- 
0107,719,. 
1665,723'i 
01C9,72'67 
1082,730, 
12 78, 733-, 
10.85,737, 
D4HT74P7 
1096,744, 
1247,747, 
1199,751, 
ITS 7, "7 54 T 
1203,758, 
1190,761, 
0579,765, 


0120,706,0117* , 
■ 1210 * 

:0 12 177 137“C 1"8 O’T— 
1231* 

,0122 , 720, C151* 
1664* 

■0v4l8 7 7277l3e2*"— ” 
1081* '4 , 

0421,734,0417* 
1220 * 

T3 43-T74'1T1‘35'6-* 

1223* 

1246 , 748,0351* 
1313* 

lies-, '75571232* 

12 80* 

1,19 1,762,109?* . 

119 7* . 


,F 120 U0,1202,'0* 
F53, .16* 

P601, 1232,602, 1$ 


, . OUTPUT ■ 

, TIME iALT,V( REU ,V(AI R) ,t. . 

V (IN ), ALPHA ( ACT) , ALPHAUUfiV) , 

“RfHAS F , R { EC ),7 G AM { R GD )“, GAM crGC') 

G'AM:( IGD) , PH I (ACT) ,PHI (COM),' 

EONG ',LAT (GD) ,AZ(RGD) ,AZ ('IG’C) 

AZll GO) ,DELE,DELBF 

-MaATTGC) ,ey,CL — 

NT,HMI6,HMD'F - 
. OB.AR ,HREFTD ,CZ,CD(CX) . 

N;T:(;E ST) ,HMQE ,'DELSB(HL) , . . 

~ 1:10 TE RR ',-R E , R'H D V ( N ) a 1 0 

0-sbJ:at,hmsb ,0FLsa (wu 

0'0€RR,HREF, HERROR,RHO 
N2C, DRAG, HOOT.'. • 

XR'W, YRW , ZRW , I WOW : 

LWRL ,IRET6F ,HRFD0T , ^ 

. XDGTRW, YDGTRW,ZDnTRW,DP 
, CPV, R(GC ) ,K (TG ) ■ 

AUTO PITCH 6 ROLL IND. 
INTEGRATION TIME STEP- 
PHASE 8 CASE STOP CONDITION 
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r 


x:ir2^»0$ 


1356 L/D 

“1 197 — RANGT^Cr-G'O-NMT 

1210 ALPHA COM ~ DEG 
100 TIME - SEC 

120 V RELATIVE FT/SEC 

~I0 8- 'ALTITUDE -~FT 

^^02 DYN PRESS LB/FT2 
166A RF DEFLECTION- DEG 
120S NOP. LOAD CQM-G 
"166 5* 'EL DEFl.ECTTON-“DEG” 
1310 ALPHA - DEG 
1231 BANK angle - DEG 
1380 BANK ANGLE - DEG 
614 '"LOAD- FACTOR G 

121 GAMMA - DEG 

IC81 HINGE MOM BF IN-LB 
1082 IN EL HM - IN-LBS 
"1. 0 6 3 — 0 U'C'-OL—HM “~1 N- LB'S' 

1187 CROSSRANGE - FT 
1186 DOWNRANGE - FT 

1188 H ABOVE RW-FT 

"1247 — HPEF -FT — 

1223 SB OEFLECTION-OEG 
1096 SB HM - IN-LBS 


>7i5r~ 


• X IOOt 
XIOQ, 
-X 100, 
X 100, 
XIOO, 
X 12 0, 
"X120, 
X120, 
X120, 
X120, 
"X120, 
X 120, 
X 120, 


' A 

>^i20, 

Ms XI 20, 




'X120, 
X120, 
X 120, 
X I2C, 
XI 186 
X 120, 
X 120, 
X 120, 
'X lO’O, 
XIOO, 
X 120, 


1356,-1 
0iC8,-2 
0120 ,> 1 - 
0402, -1 
01 . 21,-1 
0108.-1 
1665,-2“ 
16 64, + 1 
1081,-1 
1082,-1 
1083,-1“ 
1310,-1 
1210,-5-0 
1380,-1 
l2 81', + 0“ 
1223,-1 
1096,-1. 
1209, -i 
, I18 7', -r 
1188,-1 
12 47,0 
1197, -1 

ii 87 ,-r~' 

1186,0 

0614,-1$ 


•SX40 20 

X9 , 0, 800$ 


ECF 
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A-ARRAY LOCATIONS OF ARGUMENTS, RESULT AND POST MULTIPLIER 

TAPE NO. XI 4728 


DATA 

X-ARG 

Y-ARG 

Z-ARG 

RESULT 

POST-MULTIPLIER 

FILE# 

SIZE 

DESIGN HEAT- LOADS 

1655 

- 

- 

1667 

587 

17 

55 

NORMAL SHOCK VISCOSITY 

1660 

1661 

- 

1657 

588 

18 

251 

LOCAL HEAT RATIOS 

1‘655 

1600 

- 

1455 

1390 

19 

465 

TRANSITION REYNOLDS NO. 

1655 

1600 

> 

1452 

587 

20 

465 

INSULATION THICKNESS 

1655 

1658 

- 

1656 

587 

27 

575 

PAN^L AREA 

1655 

- 

- 

1652 

587 

28 

55 

correlated INV. Cd dis. 

1060 

401 

- 

1061 

1030 

60 

83 

CORRELATED INV. Cj_ DIS. 

1060 

401 


1062 

1031 

61 

83 

CORRELATED VIS. C^^ DIS. 

1060 

1353 

- 

1061 

1032 

62 

53 

CORRELATED VIS. C^^ DIS. 

1060 

1353 

- 

1062 

1033 

63 

53 

NORTH/SOUTH WIND PROFILE 

1064 

108 

- 

335 

1026 

64 

30 

EAST/WEST WIND PROFILE 

1065 

108 

- 

336 

1027 

65 

30 

DESIGN DENSITY PROFILE 

1063 

108 

- 

339 

1029 

66 

113 

ALTITUDE RATE ERROR 

1054 

120 


1055 

1034 

67 

99 

RAN6E-T0-G0 ERROR 

1056 

120 

- 

1057 

1035 

68 

87 



COEFFICIENT 


VISCOUS LONGITUDINAL 




ACm^ . 


SEGMEfnED AND REARRANGED 0V102 AERO DATA TAPE INFORMATION, 

TAPE NO. X17990 


GEMASS ARGUMENTS 
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DATA BOOK pyi p ! KfcT7tr 
REFERENCE* 


3,n-n,l2,13 1/82' 123 

8 ( 

TABLE 3.11-2 

1 ' 2/83 123 


3/84 • 123 


3.n-3a,b & 4/91 

TABLE 3.11 -3 


3.n-4a & 
TABLE 3.11-4 


7/88 

785 

8/89 

785 

9/90 

785 


NUMBER OF 
M VALUES 


MACH RANGE 
IN FILE 



NOTE: Data Book Reference stands for Figure Number unless otherwise noted 

jNOTE: Number after / represents corresponding fi]e number on basic 0V102 Aero Data Tapes 













COEFFICIENT 


AERO UNCERTAINTIES 

ACLyy 

ACpjy 

A^yy 

A^L'/o^y 

acdv„___ 


— iiu 

-^V 
A^nia . 


tiiciusure (! 
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SEGMENTED AND REARRANGED 0V102 AERO DATA TAPE INFORMATION 

. TAPE NO. XT 7990 


GEHASS ARGUMENTS] BOOK 

REFERENCE * 


X I Y I Z 


FILE NSIZE 

•ff-ff j 


3.12- 1 

3.12- 2 

3.12- 1 

' 3.12-12 

3.12- 13 
•3.12-n 

3.12- 10 

3.12- 9 

3.12- 9 

3.12- 6 


10/117 . 
n/1,18 

12/119 

13/131 

14/132 

15/133 

16/146 

y I 

17/147 

18/J48 

19/128 


•• NUMBER OF 
M VALUES 

1 

6 


MACH RANGE 
IN FILE' 


.25-10.0 


,25-10:0 


* NOTE: Data book reference stands for Figure Number unless otherwise noted 

** NOTE: Number after / represents corresponding file number on basic 0V102 Aero Data Tapes 
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COEFFICIENT 


SEGMENTED AND REARRANGED 0V102 AERO DATA TAPE INFORMATION 

limS ARGUMENTSI rZZTZZY 


X j Y ■ Z 


REFERENCE 


* FILE NSIZE 


NUMBER OF 
M VALUES 


MACH RANGE 
IN FILE 


HINGE MOMENTS 


a GgB ^ 3.6~'lB-20 


20/29 


1 . 2 - 10.0 


3.6-3D-31 


21/62 


3.6-2-17 


22/77 


I (5e • t 3.6-2-17 


23/78 . 1059 


BASIC LONGITUDINAL 


SPEEDBRAKE 

'^'"Lsb 


M a 



TABLE 3.1-1 24/1 ■ 190 


TABLE 3.2-2 25/2 190 


TABLE 3.3-1 26/3 190 


3.5-51-66 


3.5-51-66 


3.5-51-66 


27/19 748 


28/21 748 


29/23 748 


1.2-20.0 



* NOTE; Data book reference stands for Figure Number unless otherwise noted 

** NOTE; Number after / represents corresponding file number oh basic 0V102 Aero Data Tapes 
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COEFFICIENT 


BODY FLAP 
AClbp 


ELEVON EFFECTIVENESS 
ACt,, 


BASIC LONGITUDINAL 
C, 


SPEEDBRAKE 

ACi 

i-r r> 


ACd^b 

ACmcD 


SEGMENTED AND REARRANGED OVl02 AERO DATA TAPE INFORMATION 

tape" Normsoo 


GEMASS ARGUMENTS 


DATA BOOK • 
REFERENCE * 

FILE 

irk 

NSIZE 

3.5-84,86,88 

30/39 

379 . 

3.5-84,86*88 

. 31/40 

379 

3.5-83,85,87 

32/41 

379 

3. 5-1 -7 

33/12 

1377 

3. '5-8 -16 

34/13 

1377 

3.5-17-22 

35/14 

1377 

TAbUE 3.1-1 

36/1 

173 . 

TABLE 3.2-2 

37/2 

173 

TABLE 3.3-T 

38/3 

173 

3.5-51-66 

39/19 

748 

3.5-51-66 

40/21 

748 

3.5-51-66 

41/23 

748 


• NUMBER OF 
M VALUES 


* NOTE: Data book reference stands for Figure Number unless otherwise noted 

** NOTE: Number after / represents corresponding file number on basic 0V102 Aero Data'Tapes 



I 




I 


MACH RANGE 
IN FILE 


1.2-10.0 


1.2-10.0 


,25-1 .3 


.25-1.3 
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SEG^CTED AND REARRANGED 0V1Q2 AERO DATA TAPE INFORMATION 

TAPE NO. XI 7990 
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SEGMENTED AND REARRANGED 0V1Q2 AERO DATA TAPE INFORMATION 

TAPE NO. XI 7990 
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SAMPLE INPUT INSTRUCTIONS 


INSTRUCTIONS FOR CENTRAL COMPUTER COMPLEX COMPUTER RUNS 

too NOT FILL IN shaded AREASf WB5 /.2 


PROCRAMMCR SCOMMCNIS 


1 10^ EXEcIL 


QADGc NO 

sox NO 

PHONE NO DATE 

M70&^3 

SETA 


€sr TiMC^ 

MAX TIME 

PAGES OUTPUT SEG NCL 

S' 

JO 

C> 


DT-! 


OPERATING SYSTEM 

TYPE OF RUN 

I10S6XEC1I [j/j’ 

32COSCOPE -| ~j 

PBOO. | >-^ TEST 

lice EXEC VIII [ 1 

3200 SMARTS j ] 

OTHER (EXPLAIN 8ELOWI 

iiogcoaoL I 1 

32000THER 1 1 


INPUT TAPES 

yiORKI^.Q 

OtJTPUT TAPES 


NO NO. 

FASTR DRUM 

FILES files 


S' O 


psrmanent 
fASTRANO FILES 


UNIT J. RE£LNO 



MICROFILM PLOT LABEL 

JSC Form 2079 (Rev Dec 7.4) NASA-JS: 


RAhlPV ' 

i^AjBSL. 


COMMENTS 


ROOM 

DATE 

SEpI 

y/26/77 

TAPE DENSITY 

NO. COPIES 


/ 

jPROCeSSJNG REQUESTED 

^O&o 



OPT- 1 


ACTUAL 7 we U 5 Ace 


(aJbs / . 2 /^T/JcroRy 

OPERA now USE ONLY 


SEQUENCE NO. REEL NO, 


ABNORMAL STOPS 


EXCESS OUTPUT I ] 


EXCESS time: [ I 


OTHER (EXPLAIN QEUOW) [ 


wamUmKM 

































ORIGINAL PAGE IS 
OF POOR QUALITY 
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^l*OQR__QrrAirm^. 




























QUALrr y 


88 
































































































































iftiGiNAL Page t<? 

POOR QUalit^ 
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ORIGINAri page B 

OP POOT? ATT A T rrtvtjf 

















ORIGINAL PAGE IS 
OF POOR OTIATJTY 
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OEIGINAL PAGE IS 
Of POOR QUALITY! 
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OEIGINALP^^ 
^ POOR OUAUi-Xl 
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FILE ST-RUCTUPsBS 


TAPE 

GEMASS X07475 

ENTRY-TO-TOUCHDOWN 
SIMULATION 

-X06938 


FILERS ) 
1 OR 2 

1 OR 2 


AERO DATA 
TABULAR DATA 


XI 7990 1 TO 57 

XI 4728 17 TO 20, 27, 28, 

AND 60 TO 68 ■ 


-no- 



ENCLOSURE (10) 



PROGRAMMERS 


DR. F. GARCIA, JR, 

NASA/EX4 

PH. 483-5557 

K. M. NORRIS0N 

NASA/EX4 

PH. 4,83-5557 

R. H. WAIBEL 

MDTSCO/WBS 1 .2 

PH. 488-5660 X 226 


-in- 
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POTENTIAL CORE REDUCTIONS 


ROUTINE 


OCTAL DATA STORAGE 

OCTAL CODE STORAGE 

RADAR 

NOT 

424 

563 

HOUR 

REQUIRED 

31 

170 

'ORBIT 

FOR 

144 

270 

■R0CKM2 

1 

ENTRY-TO 

17 

111 

i 

BODY M2 ’ 

TOUCHDOWN 

23 

33 

CLCT . 

SIMULATION 

37 

720 

TPS 


537 

1653 

SBOOM 


100 

276 

DAPGEM 


147 

1230 

AUTLNO 


270 

2302 

TAEMG 


313 

2613 

WIND 


45 

201 

CONGIP 


306 

2541 

GIDNAV 


15 

64 

EGRT 


137 

735 

TAERO^ 


166 

1236 

-V.ECW.IN 


2303 

2442 

EXECTl 


AT LEAST 1671^ 

?3 

EXECT5 




EXECT8 


?5 
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NOTES: 
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sample of DUKMrV- ROUTINE 

FOR RADAR, RADAR 

SUBROUTINE RADAR 

RETURN 

-END 


(1) Additional 1181 decimal locations available in A-array. 

(2) See Enclosure 6 for additional A-array locations .made 
available by not using Aero Data Tape. 2966 decimal 
locations available tor not .using hinge moment 
calculations. 

(3) Overlain; reduction determined by other program mods 

in use; not required .tor-entry-to-touchdown rsimulati.on.; 
controls use of subprogram 44. 

(4) Same -as 3 except controls subprograms 21 , 22, 23, and 24, 

(5) Same as 3 except controls subprogram 35. 
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'APPENDIX I 


DISPERSIONS 

The ability to run several dispersion sources (i.e., off-nbminal 
conditions) with an entry-to~touchdov/n trajectory have. been included 
in the simulation. These are; 

T) Guidance., Navigation, and Control error 

2) Atmosphere (neutral) 

3) Winds 

4) Aerodynamic errors 

Some ,of these existed with previous forms of the GEMASS simulation, 
but more options are available with the current program. Also, a data 
tape (X17428) has been provided that contains tables of dispersion 
values. The following ‘discussion of the error sources deal with the 
options available and how to use the tape. Enclosure 6 gives- additional 
information on the use of this tape. 
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GN&C ..Error Sources 

In order to simplify the running of dispersed trajectories with 
GN&C error sources a number of modifications to GEMASS were made. 

The error sources that have been identified as attributable to 
guidance and control are: 

1. Angle of' attack 

2. Bank angle 

3. Relative velocity 

4. Altitude rate 

5. Range to go 

6. Altitude 

The angle of attack and bank angle errors are induced via the 
respective pitch and roll dynamic equations in DAP6EM. When the 
quantities resulting from these equations are integrated to give the 
actual angles, the actual steady state angle of attack and bank angle 
will differ from the commanded angles by the desired error source 
value. Imolementation of the error sources in this way allow the 
commanded angles to remain unchanged throughout the program. Both 
error source values are input in degrees-. 'The angle of attack and bank 
angle errors are input in A(1053) and A(1052), respectively. 

The relative velocity error is input in feet per second and is 
added internally to the actual relative velocity. This modified velocity 
is then used for all guidance routines. It is Input via A(1051), 
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The capabilfty to look-up an .al.t.i tude rate error is provided in 
6IDNAV. This error is added to the actual altitude rate and this 
modified rate is used for the guidance routines. The use of this 
option is controlled by A{1054). If this quantity is greater then one. 
Table 85 will look-up the altitude rate error, A(1055). A(1034) is 

reserved as a post-multiplier. The following options are available 
from the tabular data tape: 

A(1054) = 1 The WTR error profile will be called 

= 2 The ETR error profile will be called 

In a similar manner a range-to-go error can be looked-up in 
EGRT. This quantity is then added to the nominal range-to-go in 
the guidance routines. The use of this option is controlled by 
A(1056). If this quantity is greater than one, Table 86 will look- 
up the range-to-go error, A(1057). A(1035-) -i-s -reserved as a post- 
multiplier. The following options are available from the tabular 
data tape: 

A{1056) = 1 The VITR error profile will be called 

= 2 The ETR error profile will be called 

The capability to look-up a'n afT.titudie err’or ts ’provided in 
6IDMAV, This error is then added to the actual' altitude and the 
result is used for the guidance routines. The use of this option 
is controlled by A(1040). If this quantity is greater than one, 



Table 87 will look-up the altitude error, A(104T), A(1042) is reserved' 
as a post-multiplier. No tables for this are on the tabular data 
tape. 

Note that for the altitude rate and the range-to-go errors the 
user may provide his own tables if desired. The tables on the' tabular 
data tape are set iup v/i;th the y-arguments being the profile selector 
numbers (A(4054) and A(1056)) and relative velocity as the X-ar'guments. 
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Atmosphere 

6E-MASS has been modified to allow two methods of selecting the 
multipliers to nominal neutral atmospheric quantities. This option 
is controlled by A(51) with the control card multipliers determined 
from the A-array input. The quantities affected are: 

A(350) actual atmospheric pressure 
A(351) actual atmospheric density 
A(352) actual atmospheric sound speed 
A(353) actual atmospheric viscosity 
If: A(51) = 0 no table look-up, multipliers are unity - nominal 
neutral atmosphere 

A(51) = 1, the following operations will be performed: 

A(350) = A(342) * (1, + A(338)) 

A(351) = A(344) * (1 . + A{339)) 

. A{352) = A(345) * (1, + A(340)) 

A(353) = A(349) * (1. + A(341)) 

Where: A(338) to A(341) are determined by Tables 4 to 7 respectively 
and have the following post multiplier or default values 
reserved: 

A{338): A(1037) 

A{339): A(1029) 

A(340); A(1038) 

A(341): A(1039) 
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A(342), A{344), A(345) , and A(349) are the nominal neutral 
atmospheric pressure, density, speed of sound and viscosity. 
A(51) = -1, the following operations are performed: 

A(350) = A(342) * A(338) 

A{351) = A(344) * A(339) 

A(352) = A(345) * A(340) 

A(353) = A(349) * A(341) 

Four density tables compatible with the option A(51) = 1 are 
available on the tabular data tape. The density multiplier is 
determined from Table 5 and the Y-position on this table control card 
is used to determine the density profile. A{1063) is reserved for 
this option and the following profile are available: 

A(1063)<1 70° North Hot 
= 2 ETR Hot 
= 3 70° North Cold 
> 4 ETR Cold 

The X-argument of this table is altitude and the density multipliers 
are to be used with the USSA 1962. 

If the A{51) = -1 'option is chosen, a table of multipliers fn 
the form described for this option must be provided along with a 
referencing table control card. Furthermore, post-multipliers of 
unity should be assigned to A(1037), A(1038), and A(1039) if the 
pressure, sound speed, and viscosity profiles are not to be changed. 

If the nominal density (A (344)) is desired unchanged, a table control' 
card to assign a default multiplier of unity should be included in the 
input decks. 
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1‘Jinds 

The way 1n which wind velocity components are determined has not 
been changed from the previous GEMASS procedure; i-.e., table look-ups. 
However, several A-array -locations have been reserved for use as post- 
multipliers and profile indicators. Tables of design wind profiles 
as functions of altitude down to 6.5K feet .are on the tabular data tape. 
One table is used for the North or bputh wind profiles and another is 
used for the East or West wind profiles. The X-argument in each table 
is used to select the profile desired. Two A-array 's reserved for 
this selection are; 

A(1064) - East/West Wind Profile 
< 1 = East 
> 2 = West 

A(1065) - North/South 'Wind P-rofile 

1 = North 

2 = South 

Separate post-multipliers for each of the wind components have 
been reserved. They are: 

A(1026) X-wind component ■posit-mul.ttpTter 

A(1027) Y-wind component post-multiplier 

A{1028) Z-wind component post-multiplier 
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Aerodynamic Error Sources 

Calls 'to TABLE in TAERO allow the user add increments to both 
the trimmed C[_ and Cq determined by that routine. This option is 
controlled by the value in A~array 1060. 

If: A(1060) ^0.5 then no increments are added to Cj_ and Cq. 

If: A(f06D) > Oo the values are looked up in their appropriate 

tables and added to the nominal values. 

Tables of correlated aero data uncertainties are available on 
the tabular data tape. The Y-arguments are Mach' number or the 
viscous interaction parameter whichever is appropriate for the flight 
regime. The X-argument is A{1060) and is the series of aerodynamic 
uncertainty to be used. (See ref 5, sec. 9.1 .2.1 for the definition 
of these). The follov/ing post-multipl ier values have been reserved for 
use vp'tii thi^ option: 

A(1030): Inviscid Cp Uncertainty 

A{1031): Inviscid C|_ Uncertainty 
A{1032): Viscous Cp Uncertainty 
A{1033): Viscous Cp Uncertainty 
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